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ABSTRACT 

We present a spectroscopic analysis of HST/C(95 observations of three massive stars in the 
low metallicity dwarf galaxies 1C 1613 and WLM. These stars, were previously observed 
with 'VLT/X-shooter by |Tramper et al.| ( |2011[ 20141 who claimed that their mass-loss rates 
are higher than expected from theoretical predictions for the underlying metallicity. A com¬ 
parison of the FUV spectra with those of stars of similar spectral types/luminosity classes 
in the Galaxy, and the Magellanic Clouds provides a direct, model-independent check of the 
mass-loss - metallicity relation Then, a quantitative spectroscopic analysis is carried out using 
the NLTE stellar atmosphere code CMFGEN. We derive the photospheric and wind charac¬ 
teristics, benefiting from a much better sensitivity of the EUV lines to wind properties than 
Ha. Iron and CNO abundances are measured, providing an independent check of the stellar 
metallicity. The spectroscopic analysis indicates that ZjZQ = 1/5, similar to a SMC-type 
environment, and higher than usually quoted for IC 1613 and WLM. The mass-loss rates are 
smaller than the empirical ones by Tramper et al. ( 2014| l, and those predicted by the widely- 
used theoretical recipe by |Vink et al.| ( |2001| l. On the other hand, we show that the empirical, 
FUV-based, mass-loss rates are in good agreement with those derived from mass fluxes com¬ 
puted by Lucy ( 2012) l. We do not concur with Tramper et al. (2011 2014) 1 that there is a 
breakdown in the mass-loss - metallicity relation. 

Key words: stars: early-type - stars: massive - stars: mass-loss - stars: wind - stars: funda¬ 
mental parameters - galaxies: individual: IC1613, WLM 


1 INTRODUCTION 

Hot star winds are radiatively-driven by the transfer of photospheric 
photon momentum to the atmospheric material through absorption 
and scattering by spectral lines jCastor et al.||1975| jVink et al.| 
|2001| l. The basic properties of stellar winds therefore depend on 
both the number of metal lines available to absorb photon momen¬ 
tum and on their ability to absorb (i.e., their optical thickness). 
CNO and intermediate elements are dominant line drivers for the 
outer supersonic part of the winds, while iron group elements are 
responsible for the inner, subsonic part where mass-loss rate is set 

* Based on observations made with the NASA-ESA Hubble Space Tele¬ 
scope (program GO 12867), obtained at STScI, which is operated by 
AURA, Inc., under NASA contract NAS 5-26555. 

f Based on data products from observations made with ESO Telescopes at 
the La Silla Paianal Observatory under programme ID 085.D-074I 
f E-mail: Jean-Claude.Bouret@lam.fr 


I Vink et ^|200l) . The mass-loss rates, M, of massive hot stars 
are therefore expected to depend on the metal content Z. Based on 
Monte-Carlo simulations using an extensive line list, |Vink et al.| 
1 2Q0l| established a theoretical relation, M oc 
the observational side, [Mokiem et al.| P007| l derived a relation 
for early-type stars in our Galaxy and in the Magellanic Clouds 
(Z/Zq = 0.5,0.2 for LMC and SMC, respectively), obtaining a 
relation M oc that is consistent with the radiatively- 

driven wind theory. 

These conclusions have been recently questioned by |Tram- | 
Iper et al.|(|2011|(, who observed six O-type stars in galaxies with 
sub-SMC metallicities, Z/Zq = 1/7 (IC 1613, WLM, and NGC 
3109) with Ylfl/X-shooter. From the modeling of the optical spec¬ 
trum, they derived mass-loss rates that are higher than the values 
expected from the theoretical relation for the stars metallicity. The 
derived mass-loss rates are more consistent with LMC-type metal 
abundances. Several uncertainties were investigated that could af- 
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feet these determinations, including the metallicity of the galaxies 
(measured from the 0/H ratio) and the correction for nebular emis¬ 
sion, or the influence of stellar multiplicity on the derived mass- 
loss rate by dilution or wind-wind collisions. None of the solutions 
considered by |Tramper et al.H2011^ , however, could explain their 
unexpected findings. These original results have been further sup¬ 
ported by the analysis of an extended sample of 10 stars (|Tramper| 
let al.|2014| l. 

If confirmed, such a breakdown of the mass-loss rate scaling 
with metallicity would have dramatic consequences on the ultimate 
evolution of these objects, hence far more serious consequences 
than a “simple” revision of our general physical understanding of 
massive star properties. Indeed, the lower mass-loss rates that are 
predicted at low metallicities imply that the removal of angular mo¬ 
mentum during a star’s evolution is substantially reduced at lower 
metallicities, increasing the effects of rotation, especially on mix¬ 
ing jMaeder & Meyn^|2001^ . How to correctly treat the effects 
of rotation and metallicity on mass-loss and mixing is relevant to 
any field that needs to understand the life cycles and feedback of 
massive stars. These stars are the primary source of production of 
all elements heavier than O (except Fe), that they release in the in¬ 
terstellar medium via their stellar winds or when they explode as 
core-collapse supernovae. The chemical and dynamical evolution 
of galaxies is thus set by massive stars ( |Freyer et al.||2003) . An¬ 
other example is the collapsar model for long-duration gamma-ray 
bursts ( |Woosley|1993) , which requires evolved massive progenitors 
with large amounts of angular momentum and mixing, possibly to 
the point of completely mixing the star to a homogeneous compo¬ 
sition ( |Yoon & Langer|2005]|Martins et al.|2013^ . Furthermore, the 
WR/O star ratio in classical single-star evolution models is also ex¬ 
pected to be modified by the metallicity dependence of the winds: 
low-metallicity environments should produce less WR stars, which 
therefore impacts the number of type Ib/c supemovae (jMeynet &| 
|Maederf2005^ . On the contrary, should stronger winds prevail at 
low metallicity, this would imply lower numbers of collapsars, and 
of type Ib and Ic supemovae in the high-z Universe. 

|Tramper et al.| ( [2011[|2014l ) based their analysis of the mass- 
loss rate mostly on the Ha line. Several studies have shown that the 
sensitivity of Ha is limited to winds with mass-loss rates greater 
than about 10“^ Mq yr“^, and this line only shows core filli ng for 
mass-loss rates between 10~® and 10~^ Mq yr“^ (see e.g. Mar- 


Icolino et al.|2009^ . Additionally, wind clumping was neglected, al¬ 
though the presence of clumps has been established even at low 
metallicities (see |Pul s et al.|200 ^ and refs therein). [Tram per et al.| 
( |201 1[ |2014^ argue, however, that they expect that their conclu¬ 
sions would remain unchanged, should clumping be included in the 
analysis. This statement is quite surprising since, Ha being sensi¬ 
tive to the square of the wind density, neglecting clumping leads to 
overestimate the mass-loss rates by a factor 1/-^// (where / is the 
volume-filling factor associated with clumping). This corresponds 
to factors of 3 to 10 that are of similar amplitude than the expected 
metallicity effect between a Galactic massive star and a counterpart 
with the exact same photospheric parameters in the SMC typically. 
This raises the question of the reliability of the optical diagnostic 
for mass-loss determination, as well as of the underlying metallic¬ 
ity adopted for the galaxies IC 1613, WLM and NGC 3109. 

Far-UV (FUV) spectroscopy, on the other hand, is the ade¬ 
quate tool to fully address and resolve this outstanding issue of the 
dependence of hot, massive star mass-loss rates with metallicity. 
Spectroscopy in the 1150 - 1800 A waveband provides many of 
the major wind line diagnostics, including the N V AA1238, 1242 
doublet that is also temperature, abundance, and X-ray ionization 


dependent, the C IVAA1548, 1550 doublet that is particularly useful 
to diagnose weak winds, and the Balmer alpha transition of He II 
(A1640) that provides an extremely valuable clumping and abun¬ 
dance diagnostic for stars with strong winds ( |Hillier et al.|2003) ). 
In addition, N IVA1718, Si IVAA1395,1402 (also luminosity sensi¬ 
tive) and O VA1371, are also sensitive to the clumpy nature of the 
winds of early O stars. Many other useful diagnostics are also found 
in this spectral range, such as Tgff indicators from C, N, O, and Fe 
lines that can be used instead of the classical helium ionization bal¬ 
ance jHeap et al.|2006] ). Several lines are also good indicators of 
CNO and iron abundances (e.g. |Hillier et al.||2003] |Bouret et al.| 
|2003| l. The numerous iron lines (Fe III to Fe V) present in the FUV 
spectrum provide a direct measurement of the stellar metallicity. 

The goal of the work presented here was to corroborate (or 
not) the breakdown of the mass-loss rate - metallicity relation at 
sub-SMC metallicities with FUV observations with HST/C05 of 
the three brightest stars from the original study by |Tramper et al.| 
( |2011^ . The paper is organized as follows: we present the observa¬ 
tions in Sect. 1^ and we discuss the spectra morphology in Sect.|^ 
We describe the model and method we used for the determination 
of the stellar parameters in Sect. we present the results of this 
modeling in Sect. Then, in Sect. we discuss what the actual 
metallicities of the galaxies are, while a possible binary status for 
the targets is presented in Sect.|^ Systematic uncertainties in the 
mass-loss determinations are discussed in Sect. and wind prop¬ 
erties at low metallicity in Sect.|^ In Sect.[T^ we discuss how the 
wind properties we derive compare to theoretical expectations. We 
summarize our conclusions in Sect. 


2 OBSERVATIONS 


We used the Cosmic Origins Spectrograph (COS) onboard the Hub¬ 
ble Space Telescope (HST), to obtain high-resolution (R « 20,000) 
FUV spectra of the three brightest stars from |Tramper et aL] ( |201 l[ l, 
namely A13 and Bll in IC 1613 and All in WLM (see also 
|Bresolin et al.|2006l |2007| >. Aside from their intrinsic brightness, 
oth er selection criteria were tha t i) IC1613-A13 is the only star in 
the Tramper et al. (2011 2014 ' sample with a mass-loss rate (M) 


consistent with the theoretical prediction; ii) Al 1 is one of the only 
two stars in |Tramper et al.| ( [2011^ with a mass-loss rate unambigu¬ 
ously larger (~ an order of magnitude) than predicted, but the other 
object is almost one magnitude fainter; Hi) All and Bll are late- 
type O supergiants (~ 09.5 I) with much cooler effective tempera¬ 
tures (Tgff « 30, 000 K) than A13. Interestingly, M for B11 could 
be compatible with theory within the errors. On the other hand. All 
is notably brighter than Bll (0.35 dex), has lower surface gravity 
(0.2 dex) and a higher mass-loss rate (0.6 dex), suggesting that it 
is more evolved than Bll. Observing both stars in the FUV with 
HST/C05 would allow us to further constrain their evolutionary 
properties and derive unambiguous values for the mass-loss rates 
and address the issue of clumping in low metallicity stars. 

Our program GO 12867 (PI: T. Lanz) was granted 18 orbits. 
Each target was observed in a sequence of 6 science exposures with 
COS using gratings G130M (A1291, A1327) and G160M (A1577, 
A1623). This sequence provides coverage without gap of the FUV 
spectrum between 1132 and 1798 A. Overlapping segments have 
been co-added to enhance the final signal-to-noise ratio (SNR). 
Thanks to the high efficiency of COS in the FUV and little ex¬ 
tinction toward the host galaxies, we achieved typically SNR « 25 
per resolution element, in 6 HST orbits for each target (2 orbits 
for G130M and 4 orbits for G160M; the latest grating covers the 
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C IV AA1548, 1550 doublet and He II A1640). Informations about 
the targets, including UBVRl photometry and exposure times are 
summarized in Table[T] 

COS spectra are flux-calibrated (see |Massa et al.|2013[ for de¬ 
tails and limitations) providing the flux distribution in the FUV 
wavelength range. We shall use the spectral absolute fluxes, to¬ 
gether with UBVRl photometry (Table to construct spectral en¬ 
ergy distributions (SEDs) and constrain the intrinsic stellar lumi¬ 
nosity (see Sect. |4.1^ . It is possible, however, that the spectra of the 
three targets are affected by physical or line-of-sight, unresolved 
multiplicity. Indeed, the nominal primary science aperture of COS 
has a radius of 1.25 arcsecQ At the distances of IC1613 and WLM 
(cf. Table [TJ, this subtends a field of view with physical size of 
9-12 pc across, several times the size of a cluster like R136 (e.g. 
ICrowther et al.|2010t . 

On the other hand, the FUV range of a spectrum is dominated 
by the hottest objects. Shall hotter objects than the targets be in the 
COS entrance aperture, the spectra will generally present signs of 
contamination (e.g. lines from higher ionization stages) which is 
not seen (see Sect. |^. More importantly, if the stars are in binary 
systems and their mass ratios differ greatly from uniiA, then the 
derived physical properties reflect those of the primaryH 

Furthermore, ground-hased observations with spatial resolu¬ 
tion indicate that the two targets in 1C 1613 are isolated within the 
COS field of view. We note also that |Tramper et al.| p011| [2014^ 
found that their luminosities are compatible with those of (single) 
stars with same spectral type/luminosity class in the MCs (e.g. 
[Massey et al.||20()^ . Therefore, we are confident that the physi¬ 
cal parameters we can derive from the modeling of the FUV (- 1 - 
optical) spectra truly reflects the properties of the targets. The sit¬ 
uation is not as clear for WLM-Al 1, where a significant source of 
optical flux is present within a 2.5” diameter circle (hence the “a” 
subscript of the luminosity class) in optical images. Since the flux 
level in the FUV, plus the whole UBVRl photometry, is also quite 
high, it might indicate that this target is actually a binary system 
with components of similar spectral types (more exactly similar ef¬ 
fective temperatures). We will come back on this point in Sect.|^ 

We also retrieved the optical, WlT/X-shooter spectra of the 
three targets (first presented in |Tramper et al.|2011) from the Phase 
3 ESO archives. These spectra have been processed using version 
v2.2.0 of the X-shooter pipeline. We used the spectra obtained in 
the UVB (300-550 nm) and VIS (550-1020 nm) arms with slit 
width of 0.8” and 0.9”, which provide a spectral resolving power 
R = 6200 and R = 7400, respectively. Typical SNR per resolution 
element in the UVB spectra are 25, 20 and 18 for 1C 1613-A13,1C 
1613-Bll and WLM-Al 1, respectively. In the VIS arm, SNR per 
resolution element ranges from 10 to 12. 

Following [Tramper et al.| ( |2014| l, we corrected the spectra 
from the nebular emission in the B aimer lines (and He I for 1C 
1613-A13), with special care for the Ha line, as this is the prin¬ 
cipal mass-loss diagnostic available in the optical. These spectra 
have been used in the modeling process to better constrain the pho- 
tospheric parameters (especially Tgff and log g) and to check how 
the synthetic Ha profiles produced by FUV-based M compare to 
observations. 


^ Note that the COS handbook ^Holland et al.|20T^ mentions that “the 
aberrated beam entering the aperture allows objects up to 2 arcsec from the 
center of the aperture to contribute to the recorded spectrum” 

^ But see R136A for a case not showing strong evidence, from spec¬ 
troscopy and the SED, for binarity even though its a multiple system. 


3 SPECTIfA MORPHOLOGY 

In a first step, we compared the COS spectra of 1C 1613-A13, 
IC 1613-Bll and WLM-Al 1 with existing FUV spectra of SMC, 
LMC and Galactic stars of similar spectral type and luminosity 
class. In principle, this should provide a direct, qualitative, model- 
independent, indication about the mass-loss rate - metallicity re¬ 
lation. The line of reasoning is straightforward. The strength of a 
P Cygni profile is directly related to the optical depth of the wind 
in the corresponding line. The optical depth is related to the total 
number of absorbing ions in the transition giving rise to the spec¬ 
tral line and can be written as Trad oc ALqiAe with Ae and gi the 
abundance of element E and its ionization fraction for stage i. It is 
therefore readily seen that, for stars with similar physical proper¬ 
ties and evolutionary stages (as should be for stars within the same 
spectral/luminosity classes), hence yielding similar ionization frac¬ 
tions, Trad scales with the mass-loss rate and the abundance (which 
is expected to scale with metallicity). 

We selected stars with spectral classification as close as pos¬ 
sible to those of IC 1613-A13, IC 1613-Bll and WLM-Al 1 al¬ 
though a perfect match was not always possible. Also, whenever 
possible, we chose stars with high resolution spectra, such that 
we can further compare the numerous lines from iron ions. High- 
resolution lUE spectra exist for Galactic massive stars of all spec¬ 
tral type and luminosity class. In the Magellanic Clouds (MCs), 
several lUE spectra are available but they have poor signal-to- 
noise ratios and were not considered for this comparison. In the 
MCs, good signal-to-noise high-resolution spectra have been ob¬ 
tained with HST but for a limited number of objects. For early- 
type dwarfs in the LMC, such good signal-to-noise high-resolution 
(Goddard High-Resolution Spectrograph, GHRS) spectra exist but 
they are restricted to stars in R136. These O-type stars often seem 
quite extreme objects in terms of stellar and wind parameters 
[Koter et aL||1998HPrinja & Crowther|1998[ [Massey et al.||2004^ ,' 

in excess of those found for Galactic stars of similar spectral type 
( [Bouret et al.|2005[|Martins et al.|2012) . We chose to use the spec¬ 
trum of a field star instead, namely LH81W28-5, whose spectrum 
was obtained at low resolution though (G140L) with STIS ( |Massey| 
|et al.|20d4l ). We could not find a high-resolution spectrum ( COS, or 
STIS or GHRS) for LMC late-type supergiant and we had to rely on 
the more modest resolution provided by EOS for BH70; we use it 
for comparison to both IC1613-B11 and WLM-Al 1. We could find 
only one late-type supergiant in the SMC, namely AzV 327, with 
high-resolution (STIS) FUV spectrum. Here again, this star is used 
for comparison to Bll and All. The IDs of the stars we used for 
the morphology comparisons are summarized in Table together 
with additional indications about their FUV spectra, V-magnitude 
and reddening. 

For the comparisons presented in Fig. [T] (Bll and All) and 
Fig.[^(A13), we used absolute fluxes, corrected for reddening, such 
that we avoid potential problems related to the choice of a contin¬ 
uum to rectify the spectra. Fluxes for Galactic, LMC and SMC stars 
have been scaled by the appropriate ratio to match the continua of 
the three targets in spectral regions of interest (i.e., around wind 
profiles). These scaling factors present relatively small deviation 
from those that would be obtained from a scaling to the distance 
modulus of IC 1613 and WLM as given in |Tramper et al.|(2011| |. 

We first discuss the cases of IC 1613-Bll and WLM-All as 
they have similar spectral type/luminosity classes. Qualitatively, 
their behavior is easy to interpret in the framework of the line- 
driven wind theory. The top three panels in Fig. show that the 
wind terminal velocity decreases with decreasing metal content. 
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Table 1. Basic parameters of our targets and exposure times of the COS spectra. 


Star 

Sp. Type 

U 

B 

V 

R 

I 

E{B - V) 

dist. 

(kpc) 

Mv 

Texp (s) 

G130M G160M 

IC 1613-A13 

03-04 V((f)) 

17.66 

18.73 

18.96 


19.26 

0.03 

721 

-5.42 

4640. 

10006. 

IC 1613-Bll 

09.5 I 

17.45 

18.49 

18.62 


18.78 

0.08 

721 

-5.92 

4581. 

9927. 

WLM-Al 1 

09.7 la 

17.23 

18.27 

18.38 

18.45 

18.57 

0.04 

995 

-6.73 

4625. 

9992. 


Note : Spectral types ai‘e from Tramper et al.|(2QTT) and |Garcia et al.|j2014^ , while the photometry is taken from |Garcia et al.|j2009||2014^ for the stars in 
IC1613 and from|Massey et al.| 2007| for WLM-Al 1. Color excess E{B — V) is calculated using the calibration byjMartins & Plez j2006| for Galactic Stars. 
Absolute magnitude are calculated using the quoted distances to IC1613 and WLM from |Pietrzyhski et al.H20Q^ and |Urbaneja et a .^2008) respectively. 


IC1613-B1 1 



Wavelength (A) 


WLM-Al 1 



Wavelength (A) 


Figure 1. P Cygni profiles for Bll (left) and All (right) in black. Overplotted in red, green and blue are the spectra of the Galactic, LMC and SMC stars 
(respectively) chosen for comparison (cf. Tablej^. For plotting puipose, the fluxes were scaled by a factor 10^'*. 


from Galactic to IC 1613AVLM, while the P Cygni profiles get si¬ 
multaneously weaker. This especially clear for the second and third 
panels for each star, which demonstrate that this trend holds for 
the resonance doublets of Si IV and C IV. On the other hand, cau¬ 
tion should be used concerning the N V resonance doublet at 1238, 
1242 A, which is notoriously sensitive to X-ray irradiation, that is 
not known for the stars used for this comparison. P Cygni profiles 
are well-defined for the Galactic and LMC stars but only a weak (if 
any) emission component is observed for the SMC star, as well as 
IC1613-B11 and WLM-Al 1. The C IV AA1548, 1550 profiles still 
show a blueward extension that is related to the wind whereas it is 
barely seen on the Si IV doublet because of the ten times lower 
abundance of Silicon versus Carbon. The panels showing N IV 
A1718 on the other hand, show very similar profiles for IC1613- 
B11 and WLM-Al 1 and their SMC counterpart, both weaker nev¬ 
ertheless than for the Galactic star of comparison. Overall, the mor¬ 
phology of the profiles from Fig. seem to contradict the conclu¬ 
sion by[^amperet^([20TT]lthatIC1613-Bll andWLM-All have 
mass-loss rates higher than expected for their metallicity. This of 
course assumes that the line-driven wind theory holds for stars in 
the Galaxy and down to SMC metallicity, as checked by|Mokiem| 
[eraL| ( |205^ . 

It is also easily seen on Fig.[T]that the iron lines present in the 
FUV spectra of IC I6I3-BII and WLM-AII (see e.g. the bottom 
two panels, Fe IV lines between AA1530-1570 and AA1710-1730) 
are moderately weaker, if at all, than those of the SMC star we 
chose for comparison, despite a slightly higher projected rotational 


velocity {v sin i) for AzV 327 (100 km s“^, Bouret et al. in prepa¬ 
ration). Although it could be an indication that the iron abundance 
in our targets is slightly smaller than the fiducial 0.2 Fcq of SMC 
O-type stars, it could as well be related to a lower micro-turbulent 
velocity in their photospheres. Settling this issue requires quantita¬ 
tive modeling, which we will discuss later in this paper. 

Let us now discuss IC1613-A13. As already said the N V res¬ 
onance doublet (top panel of Fig. is sensitive to the X-ray flux 
(not known for most of the stars here); it also often heavily blends 
with interstellar Lya on the blue side of the profile, which is an 
additional issue for early-type O dwarfs as their wind terminal ve¬ 
locities are expectedly higher than for the late-type supergiants. De¬ 
spite this, Voo does appear to decrease with metallicity, a trend that 
is further confirmed on the other panels. Using the C IV AA1548, 
1550 resonance doublet (third panel,), we see that Voo for the LMC 
star (« 2700 kms“^), intermediate between Voo for the Galactic 
and SMC star. Remarkably, Voo for IC1613-A13 and the counter¬ 
parts in the SMC (MPG 324) are very similar, which might indicate 
that the metallicities in both galaxies are alike. 

The conclusions are not as obvious for the mass-loss rate. For 
instance, th N V A1240 P Cygni profile in A13 is clearly as strong 
(both the absorption and emission components) as the profile in the 
Galactic star (HD 93205) and in any case stronger than in MPG 
324 in the SMC. The second panel, showing O V A1371 seems to 
indicate that, the LMC star left aside, we do observe a trend for 
lower metallicity stars to have lower M. The O V A1371 profile 
for IC1613-A13 is indeed almost completely photospheric while 
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Table 2. IDs, spectra and photometry for the comparison stars. 


Star 

Comparison Stars - Sp. Type 

Galaxy 

FUV spectrum 

V 

B 

EiB-V)*" 

IC1613-A13 

HD 93205 - 03.5 V((f)) 

IC1613 

MW 

COS GO 12867 
WE SWP07959 

7.76 

7.76 

0.28 


LH81W28-5 - 04V((f+)) 

LMC 

STIS GO 8633 

13.92 

13.74 

0.15 


MPG 324 - 04 V((f)) 

SMC 

STIS GO 7437 

14.07 

13.79 

0.07 

IC1613-B11 

HD 76968 - 09.5 Ib 

IC1613 

MW 

COS GO 12867 
WE SWP20621 

7.09 

7.21 

0.38 


BI170-O9.7I 

LMC 

EOS GO 5444 

13.06 

13.26 

0.07 


AzV 327 - 09.5II-Ibw 

SMC 

STIS GO 7437 

13.25 

13.47 

0.04 

WLM-Al 1 

HD 152424 - 09.7 la 

WLM 

MW 

COS GO 12867 
WE SWP15021 

6.31 

6.71 

0.66 


BI170-O9.7I 

LMC 

EOS GO 5444 

13.06 

13.26 

0.07 


AzV 327 - 09.5II-Ibw 

SMC 

STIS GO 7437 

13.25 

13.47 

0.04 


Calculated from B-V (this table) and (B-V)o from ^Martins & Plez|2006l . 


IC1615-A13 



Wovelength (A) 


Figure 2. P Cygni profiles of A13 (in black). Overplotted in red, green 
and blue are the spectra of the Galactic, LMC and SMC stars (respectively) 
chosen for comparison (cf. Table|^. For plotting purpose, the fluxes were 
scaled by a factor 10^^. 

the profile for MPG 324 (SMC) and HD 93250 (MW) present a 
blueward extension (stronger for the Galactic star). However, O V 
A1371 is one of the FU V line that is expected to be most sensitive to 
wind clumping ( [Bouret et al.|2005| >, together with N IV A1718.The 
later is shown in the bottom panel and the same comments hold for 
both lines (this wavelength range is not covered by the 5775 spec¬ 
trum of the LMC star). Care should be used when interpreting these 
profiles as their strengths also reflect the degree of wind clumping, 
that most likely varies from star to star. The gap in metallicity be¬ 
tween our Galaxy and the SMC is large enough however, such that 
the weakening of the wind profile caused by the factor of five in 
metallicity is visible in the profiles in this panel. Here again, the 
resolution of the LMC spectrum is too low to reach a firm conclu¬ 
sion about the strength of the O V A1371 profile relative to the SMC 
counterpart and IC1613-A13 accordingly. 

The third panel, showing the C IV AA1548, 1550 resonance 
doublet, further adds to the confusion. The expected qualitative 
trend is observed for the Galactic, LMC and SMC stars for the 
mass-loss rate (and terminal velocity). The absorption components 
of the P Cygni profiles for the three stars go from clearly satu¬ 


rated (Galactic, then LMC) to unsaturated (SMC); the emission de¬ 
creases in intensity as well. This is fully consistent with theoretical 
expectations from the radiatively-driven winds theory, but the be¬ 
havior of IC1613-A13 is at odds with this predicted trend. Although 
its Vao seems quite similar to that of MPG 324 in the SMC, the in¬ 
tensity of its P Cygni profile is much stronger, suggesting that the 
mass-loss rate is higher. One should remember, however, that the 
intensity of the P Cygni profile is determined by Mqq jy^c where 
7C IV ionization fraction of C IV and Ac the abundance of 

Carbon. Therefore, we cannot rule out that either the Carbon abun¬ 
dance or the C®"*" ionization fraction, or a combination of both, is 
really lower in the SMC star such that the P Cygni profile is weaker. 
It is also noteworthy that the C IV P Cygni profile of IC1613-A13 
shows well-defined narrow absorption components (NACs) that are 
useful to measure the wind terminal velocity ( [Prinja & Crowth^ 
|1998| l. DACs are also observed in several SMC O dwarfs of early 
spectral type ( |Bouret et al.|2013} , where the absorption component 
is never saturated (see Fig.|^ the case of MPG 324) . 

To summarize this section, we find that the P Cygni profiles 
of IC 16I3-B11 and WLM-AII are compatible with the winds of 
these stars being weaker than winds of stars with similar spectral 
types and higher metallicity. Since [Mokiem et ^ ( |2007] l showed 
that the wind quantities of massive stars in the Milky Way, the 
LMC and SMC follow a trend in quantitative agreement with the 
predictions of the radiatively-driven wind theory, it is tempting to 
conclude that the spectra of our targets also comply quantitatively 
with the theoretical expectations. The case of IC I6I3-AI3 is dif¬ 
ferent and only ambiguous conclusions about the wind quantities 
and their dependance to metallicity is obtained from a qualitative 
comparison. We will now discuss the quantitative aspects of this 
conclusion, from the modeling perspective. 


4 MODELS AND METHOD 

For a more quantitative view of the wind properties, we compared 
the observed FUV -l- optical spectra of IC I613-AI3, IC 16I3-B11 
and WLM-AI1 to synthetic spectra calculated with the NLTE stel¬ 
lar atmospheres code CMFGEN ( [Hillier & Miller||T998| ). CMF- 
GEN computes NLTE line-blanketed model atmospheres, solving 
the radiative transfer and statistical equilibrium equations in the 
comoving frame of the fluid in a spherically-symmetric outflow. 
The models we computed account for the presence of the following 
species: H, He, C, N, O, Ne, Mg, Si, P, S, Cl, Ar, Ca, Ee, Ni, repre- 
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senting more than 7500 full levels (more than 2000 super-levels) 
and ~ 135,000 lines. The density structure is described as an hy¬ 
drostatic part, in the quasi-static photospheric layers, smoothly con¬ 
nected to a /3 velocity law in the wind. The radiative acceleration 
is calculated from the solution of the level populations, and is used 
to compute a new inner structure (connected to the same /3 velocity 
law). The mass-loss rate, density, and velocity are related via the 
continuity equation. After convergence, a formal solution of the 
radiative transfer equation was computed in the observers frame 
( [Busche & Hillier|2005| ), thus providing the synthetic spectrum for 
a comparison to observations. 

We started from the stellar and wind parameters provided in 
[Tramper et al.| ( [^ll[ ) and tuned them as needed to get the best fit 
of the FUV and optical spectra together. The influence of clump¬ 
ing on the wind profiles was considered. Clumping is described in 
the framework of the optically-thin clumps formalism implemented 
in CMFGEN. The clumped wind density, p{r), is related to the 
homogeneous (unclumped) wind density p, and the volume filling 
factor / by p(r) = pjf. The filling factor decreases exponentially 
with increasing radius (or, equivalently, with increasing velocity): 
/ = /oo -I- (1 — /oo) exp(—n/uci) where Ud is the velocity at 
which clumping starts. More physical description of clumping ex¬ 
ist nowadays and will be discussed in Sect. |8.2| A limitation of this 
parameterization is that / is monotonic whereas theoretical simu- 
lations by [Runacres & Owocki| (|2002|) and observational studies 
( |Puls et al.|2006 Najarro et al.|20liy suggest it is non-monotonic. 
We found no need to tune Vd to improve the fit to observed lines , 
hence we adopted a canonical Vd = 30 km s“^ jBouret et al.|2005] >. 

The major consequence of this formalism is that mass-loss 
rates derived with clumped models are reduced by a factor \/(/oo), 
where foo is the volume filling-factor associated with clumping, 
compared to mass-loss rates obtained from homogeneous models 
(e.g. |Bouret et al.|2005] >. In the rest of the paper, we will refer to 
the “clumped” mass-loss rate (M and /oo derived from the model¬ 
ing of the FUV spectra) to compare to theoretical quantities and/or 
mass-loss rates from |Tramper et aL| ( |2014^ . 

We also accounted for the effects of shock-generated X-ray 
emission in the model atmospheres. Since no measurements are 
available for IC 1613-A13, IC 1613-Bll andWLM-All in the X- 
ray wavelength range, we adopted the standard Galactic luminosity 
ratios logLx/Lboi ~ — 7 ( [Sana et al.[|2006^ , with the bolometric 
luminosity determined in the present study. We acknowledge how¬ 
ever that there is no guarantee that this scaling ratio should pre¬ 
vail at low metallicity. |Owocki et al.| ( [2013) suggested for instance 
that this ratio should change with the physical nature of the shock 
generating the X-rays (radiative or adiabatic) and the wind density. 
Depending on the actual relation between the metallicity and the 
mass-loss rate, logLx/Lboi could change accordingly. 


4.1 Stellar and wind parameters 

The photospheric (including surface abundances) and wind param¬ 
eters have been derived following the same methodology as in pre¬ 
vious papers (see e.g. [Bouret et al.|20T^[20T3) . 

• Stellar luminosity: The distance to IC 1613 and WLM were 
taken from jPietrzynski et al.j ( [2006^ and jUrbaneJa et al.j ( |2008| l, 
respectively. We used the flux-calibrated COS spectrum, providing 
the flux distribution in the FUV wavelength range, together with 
optical photometry (cf. Table|^ to build spectral energy distribution 
(SED). The intrinsic stellar luminosity was then derived, by fitting 


the SED. Note that we used the color excess quoted in Ta ble to 
take reddening into account, using the data by jCardelli et al.|fl989| l. 

• Effective temperature and surface gravities: Tgjy was con¬ 
strained from the helium (primarily He I A4471 to He II A4542 
equivalent width ratio, as well as by fitting the line profiles) in the 
optical spectra. We add that the relative strength of helium lines 
was also used to constrain the helium abundance. Ratios of lines in 
the EUV spectra jHeap et al.|2006[ [Bouret et al.|2013| e.g.) were 
also used for further consistency, e.g. C IV A1169 to C III A1176, 
or N III AAl 183-1185 and N ill AA1748-1752 to N IV A1718 ra¬ 
tios. This is particularly important for the hottest star of the sample, 
namely A13, where the He I lines become weak (see also[Rivero[ 
[Gonzalez et al.|2012[ for the use of nitrogen optical lines to con¬ 
strain Tgjy for the earliest 0-stars). The C IV Al 169 to C III Al 176 
ratio is especially useful in early O-type stars because of its very 
good sensitivity to Tgjy and weak dependance on gravity; hence we 
used it as the primary Tgjy diagnostic in the UV. 

The Stark-broadened wings of Balmer lines, especially H 7 , were 
the primary constraint on log g. The S/N of the X-shooter spectra is 
such that it is not possible to achieve a typical accuracy (1 ct) better 
about 0.2 dex in log g. 

• Microturbulence: A constant microturbulent velocity of 15 
kms“^ was adopted for calculating the atmospheric structure. To 
calculate the emerging spectrum, we used a radially dependent tur¬ 
bulent velocity (see also [Hillier & Miller|1998f . The microturbu¬ 
lence in the photosphere was set mostly using the S V A1502 
line, that is very sensitive to this parameter over a large range of 
^eff- We further tuned it if necessary using the known sensitivity 
of some He lines to ( [Villamariz & Herrero|2000l [Hillier et al.[ 
[2003[ l. In the outer region of the wind, the turbulence was adopted 
to fit the shape and slope of the blue side of the absorption compo¬ 
nent of the C IV 1548-1551 P Cygni profile. 

• Surface abuudauce : The thousands of photospheric lines of 
iron ions (Ee III to Ee VI) present in the EUV are most useful to 
constrain the metal content of the stars. This shall provide an in¬ 
dependent check of the stellar metallicity so far inferred from the 
ISM 0/H ratio in the host galaxies of our targets. While the lines 
of Ee ions grow stronger with increased gravity and higher iron 
abundance and/or microturbulent velocity, the ratio of line strengths 
is most sensitive to temperature (see the discussion in [Heap et al.[ 
[2006[ >. Practically, we used Tgg, log g, and determined as we just 
presented to derive an iron abundance as reliable as possible. Abun¬ 
dances of a-elements included in the models (Ne, Mg, Si, S, Ni) 
were scaled down from the solar values of [Asplund et al.[ ( [2005[ l, 
adopting the same factor as for iron. Note that for the specific case 
of Si, we tuned the silicon abundance in the two late-type O super¬ 
giants, using several lines of Si III (triplet at AAA1294.54, 1298,89, 
1298.96 or the singlets at A1312.59 and A1417.24) present in the 
COS spectra. 

- Carbon abundance'. ClvA1169, ClIlA1176 are prime indi¬ 
cators for carbon abundances. These lines proved to be sensitive 
to Tgff and and C IIIA1176 further showed sensitivity to the 
wind density for the two stars with later spectral types. In this 
case, we gave more weight to C IIIA1247, which is clearly detected 
in IC1613-B11 and WLM-All and not affected by the adjacent, 
weak, N V resonance doublet. 

The WXIX-shooter spectra turned out to be quite noisy and with 
a spectral resolution somewhat too low, such that C III AA4068, 
70 were used only as secondary indicators for C/H determination 
for the late-type supergiants (Bll and All). Besides, we did not 
include C III AAA4647, 50, 51 although is clearly detected in the 
spectra of Bll and All. [Martins & Hillier[([2012t indeed showed 
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that it is very sensitive to the adopted atomic physics and other 
modeling assumption^ 

- Nitrogen abundance: In the UV, the photospheric lines N III 
AA1183-1185 and N III AA1748-1752 are used as primary diag¬ 
nostics. N IV A1718 is also used for the late supergiants, but not 
for A13, as it is clearly affected by the stellar wind in this early- 
type star. In the optical N III AA4634-42 and N IV A4058 are the 
strongest nitrogen lines in the spectrum of A13, and were used 
as additional diagnostics for nitrogen abundance. Their formation 
processes, including their sensitivity to the background metallicity, 
nitrogen abundance, and wind-strength, have been discussed exten¬ 
sively by |Rivero Gonzalez et al.H2011||2012^ . 

- Oxygen abundance: The O IVAA1338-1343 lines usually serve 
as a diagnostic for 0/H in the UV. Here however, they show a blue 
asymmetry that indicates a contribution from the wind for the earli¬ 
est star (A13) with the stronger wind, and they are very weak (if de¬ 
tected at all) in the late-type supergiants. O IIIAA1150, 1154 could 
be used instead but these lines are on the shortest side of the UV 
spectra and the flux level there is more uncertain. Similarly, O V 
A1371 is detected for the hottest star (A13) but it is very sensitive 
to the to the mass-loss rate and clumping parameters. Because the 
number of useful/reliable oxygen lines is so limited, the oxygen 
abundance was first set such that the sum of CNO is conserved, 
with respect to an initial mixture CNO scaled to the metallicity of 
the star. It was then tuned if necessary to improve the fit quality to 
OIVAA1338, 1343. 

• Wind parameters: The wind terminal velocities, I4o, were 
estimated from the blueward extension of the absorption compo¬ 
nent of FUV P Cygni profiles. Although the P Cygni profiles of 
the three targets are not saturated, narrow absorption components 
are visible, especially in N V and C IV profiles and were used to 
derive estimates of Voo, following the approach discussed in |Prinja| 
|& Crowther| ( |1998) l. The typical uncertainty in our determination 
of Voo is 100 kms“^ (depending on the maximum microturbulent 
velocity we adopted). 

Mass-loss rates were derived from the analysis of UV P Cygni 
profiles, such as C IV AA1548, 1551, O V A1371 and N IV A1718 
for A13 and Si IV AA1394, 1403 and C IV AA1548, 1551 for Bll 
and All, respectively. We emphasize again the strong dependance 
of the N V resonance doublet on the adopted X-ray luminosity, es¬ 
pecially for the stars with cooler Tgff. As logLx/Lboi is not known 
for our targets, we did not use N V A1238, 1242 P Cygni profiles to 
constrain M. 

The P exponent of the wind velocity law was derived from fit¬ 
ting the shape of the P Cygni profile. The clumping parameters, fao 
and Vci were derived in the UV domain. We used O V A1371 and 
N IV A1718 ^Bouret et al.|2005t , as well as Si IV AA1394, 1403 and 
C IV AA1548, 155. All these lines indeed present significant sensi¬ 
tivity to clumping in the temperature regime of our targets. Some 
photospheric lines in the optical presented some sensitivity to the 
adopted filling factor (and scaled M). For photospheric H I and He 
lines, for instance, this is essentially caused by a wind contribution 
that is weaker in clumped models, which produces deeper absorp¬ 
tion than smooth-wind models. 

Synthetic spectra have been convolved with gaussian instrumen¬ 
tal profiles, at the resolution delivered by HST/COS and VLT/X- 


® We note however that the influence of Fe lines found by|Martins & Hillier| 
j2012) should be less than for galactic O stars because the iron (and more 
generally metal) abundance is expected to be down by a factor of five or 
more 


shooter, respectively. Then, we further convolved these spectra 
with a rotational profile to account for the projected rotational ve¬ 
locity of the star. For simplicity, we adopted the values listed |Tram-| 
|per et al.] ( |201 l| l and checked on several FUV line profile for con¬ 
sistency (see e.g. |Bouret et al.||2013) . Given the moderate spec¬ 
tral resolving power and signal-to-noise ratio of the NYIlX-shooter 
spectra, we did not try to derive additional broadening mechanisms 
using the classical Fourier Transform method (|Sim6n-Diaz|2010). 


5 RESULTS OF THE MODELING 

The best-fit of the FUV and optical spectra of the targets are pre¬ 
sented in the Appendix, while the parameters of the corresponding 
model are listed in Table[^ We emphasize that the values we quote 
for M in this table were obtained by fitting the FUV spectra we 
presented hereabove, but they do provide consistent fits to the op¬ 
tical spectra. For comparison to parameters derived from optical 
spectrum, we used the most recent work by |Tramper et al.H2014[ l 
as the fitting procedure presented there was updated compared to 
|Tramper et al.|(|2011f to account for the scaling of Voo with metal¬ 
licity, and the qualitative conclusions already reached in |Tramp^ 
|et al.|([ 20 TT) are unchanged. 


5.1 IC 1613-A13 


We find A13 to be cooler than in |Tramper et al.|p014| >, although 
both Tgj-f overlap within the respective error bars. Since we used the 
ionization balance of C IV and C III as well as those of the iron ions 
present in the FUV to constrain the parameters, our Tgjy estimate is 
therefore less prone to uncertainties related to the weakness of He I 
lines in the optical for such an early-type O star. We note that our 
value is very close to Tgjy derived by |Garcia et al.| ( (2013[ l from the 
modeling with CMFGEN of a lower resolution COS spectrum of 
this star. The luminosity we derive, based on SED fitting, is a mere 
0.09 dex smaller than in |Tramper et al.| j2014^ . On the other hand, 
we favor a surface gravity log ^ = 3.75 over the lower value (log g = 
3.65 in |Tramper et al.|2014^ , although the error bars are quite large 
in both case. This higher value stems mostly from the influence of 
log g on the global ionization balance (see |Heap et al.|2006| ) and 
is also more compatible with the spectral/luminosity class of IC 
1613-A13. As a consequence, the stellar, spectroscopic, mass we 
derive is in good agreement with the mass quoted by|Tramper et al.| 
( |20T4t . 


The helium fraction (by number) we derive is a factor of 2.5 
times smaller than in |Tramper et al.| j2014^ . We find that such a 
high helium fraction as theirs would produce He II A1640 with a 
well-developed P Cygni profile, as well as strong wings in emission 
for the He II A4686 line (see Eig.[^. None of those are observed 
and we discarded this solution; it is also not consistent with the 
evolutionary status of this star. 

The iron abundance, which we shall now refer to as the 


global metallicity, was first set at Fe/Fe 0 = 0.14 i Tramper et al. 


|201 1| and refs, therein). Models with an iron content typical of 
SMC (Fe/Feg = 0.20) and EMC (Fe/Eeg = 0.50), as well as 
Fe/Ee 0 = 0.10, were computed. The higher and lower metallicity 
options are easily discarded but on the other hand, it is not possi¬ 
ble to rule out a SMC-like metallicity for IC 1613-A13. As already 
outlined in Sect. the strength of the many Fe III to Fe VI lines 
(as well as S V A1502) depends on the adopted abundance as much 
as on the photospheric parameters Tgjy, log g and the microturbu¬ 
lence set in the spectrum calculation (cf. Fig. [^. Within error 
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Figure 3. In black, line profiles for different transitions of iron ions (Fe Vito 
Fe IV from top to bottom) in the FUV range, observed for IC1613-A13. 
Overplotted in red, green, blue and pink are profiles obtained for models 
with the same parameters but the iron abundance, ranging to 1/2, 1/5, 1/7 
and 1/10 Fe 0 respectively. 


Table 3. Stellar parameters derived from the spectroscopic analysis 


Properties 

A13 

Bll 

All 

Spectral type 

03 V((f) 

09.51 

09.7 la 

Tgff(K) 

42500. 

30000. 

29000. 

logg (cgs) 

3.75 

3.25 

3.25 

log 7 ^ 

5.62 

5.45 

5.69 

V/gpec [M0] 

27.6ibl3. 

22.3ib7. 

53.2ib9. 

[kms“l] 

25 

15 

10 

Vtfe 

0.10 

0.10 

0.11 

ec 

7.00ib0.24 

7.00ib0.30 

7.00ib0.30 

cn 

7.78ib0.23 

7.70ib0.20 

7.70ib0.20 

eo 

7.90ib0.20 

7.96ib0.30 

7.96ib0.30 

Fe/Fe 0 

0.2 

0.2 

0.2 

[C/H] (dex)** 

-1.39 

-1.39 

-1.39 

[N/H] (dex)^ 

0.00 

-0.08 

-0.08 

[O/H] (dex)^ 

-0.76 

-0.70 

-0.70 

foo 

0.03 

O.IO 

0.10 

Vao 

2180. 

1300. 

1400. 

0 

0.8 

1.0 

1.0 

logiMpuv) 

-6.60 

-7.82'" 

-7.95*" 

^Og{Ad j'ramper) 

-5.85 

-6.25 

-5.50 

\og(Mvink) 

-6.07 

-6.51 

-6.41 

log(Miugy) 

-6.82 

-7.65 

-7.54 


Uncertainties on Tgpf are it 1000 K, it 0.2 dex for logg, it 0.1 dex on 
log L / Lq . As for the wind quantities, an uncertainty of it 0.2 dex was esti¬ 
mated for the mass-loss rates (given in Mq yr“ ^), while Voo was measured 
within it 100 kms“^. For CNO abundances, ex = 12 -I-log (X/H). 

^ Defined as [X/H] = log (X/H) - log (X/H )0 . Solar values adopted from 
[Asplund et al.|j200^ . 

See Sect. l8.3| for a discussion about the true mass-loss rates from FUV 


bars of ±5 km for these lines are well fitted either adopting 
Fe/Fe 0 = 1/5 or Fe/Fe 0 = 1/7 for the iron content. The same 
conclusion holds for the Si IV resonance doublet and the S V A1502 
line, which are purely photospheric in IC1613-A13. 

An important conclusion of our study is that IC1613-A13 
is carbon/oxygen-depleted and nitrogen-rich, compared to scaled- 
solar CNO abundances. The carbon abundance is very similar to 
the one we derived for MG 324, the SMC 04 V star we used for 
comparison in Fig. (see also |Bouret et al.||2013) l, but the nitro¬ 
gen abundance is higher. The rotation rate of A13 is only mildly 
higher than in MPG 324 (we adopted v sin i = 98 km s“^ for A13 
as in Tramper et al. 2011, while usin i = 70 kms“^ for MPG 324). 
Furthermore, we note that IC 1613-A13 is more luminous and has 
lower surface gravity than MPG 324 (and other 04 V stars in the 
SMC for that matter), which points to a more advanced evolution¬ 
ary status. An interpretation of the different abundances of these 
stars would require comparisons to predictions from evolutionary 
models, which is beyond the scope of this paper. In any case, the 
similar carbon abundance of IC I6I3-AI3 and MPG 324 implies 
that the observed difference in the C IV resonance profiles of these 
stars truly reflects a difference in mass-loss rates and wind clump¬ 
ing (Tgff is roughly the same for the two stars, hence the ioniza¬ 
tion fraction of C IV). The oxygen abundance conserving the sum 
of CNO, with respect to an initial CNO mixture that has the same 
scaling from the solar CNO mixture than the iron content of the 
model, produces a good fit to the FUV oxygen lines (O IV AA1338, 
1342 and O V A1371) and was adopted for the rest of this work. 

The major modification with respect to|Tramper et al.|j2014^ 


is about the wind parameters, that is the mass-loss rates, wind ter¬ 
minal velocities and clumping filling-factor, which we derive from 
the modeling of FUV P Cygni profiles. Fitting the COS spectrum of 
IC1613-A13 required a volume filling factor of /oo = 0.03. Aside 
from the well-known clumping diagnostics such as the O V A1376 
and N IV A1718 lines, clumping turned out to be an important 
ingredient to fit the C IV AA1548, 1550 P Cygni profile. This is 
because the ionization balance of carbon is shifted toward C^”*" 
throughout the wind, at the effective temperature of IC1613-A13 
and clumping helps C^”*" recombine to C®'*'. More generally, this is 
true for Tgff in the range usually covered by 04 dwarfs and earlier 
(typically Tpff 42,000 K). For Tgjy as high as quoted in |Tram-| 


per et al. j2014|, the ionization fraction of C^"^ is actually close to 


unity and the resonance doublet is weak for an homogeneous wind. 
Significant level of clumping is required to produce enough C®”*" 
such that a strong C IV AA1548, 1550 P Cygni profile forms again. 
This dependance of the C IV profile on the clumping volume filling 
factor might explain why its more difficult to see the trend of M 
with metallicity in A13 than in the late-type supergiants (cf. Sect 


0 . 

The clumped mass-loss rate is 2.5 xlO ^ while 

the wind acceleration is /3 = 0.8 and the wind terminal velocity is 
2180 kms“^. The later are in excellent agreement with the values 
derived by [Garcia et al.| j2014^ from the modeling of low resolu¬ 
tion COS spectrum of this star, with the SEI method (Sobolev with 
Exact Integration, |Hamann|1981[[Lamers et al.|1987[l. 
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5.2 IC 1613-Bll and WLM-All 

The photospheric parameters derived for IC 1613-Blland WLM- 
All from the analysis of their FUV optical spectra are consistent 
with their spectral classification, and in relatively good agreement 
with those of |Tramper et al.| ( |2014| l. The surface gravity we de¬ 
rive for IC 1613-B11 is actually more consistent with log g derived 
by |Bresolin et al.H2007| l. We note that the Si III lines (AA1295 - 
1299) get weaker as logg decreases. They argue against logg ^ 
3.2, which we adopted as a lower limit as it provides the best fit to 
the FUV spectrum while still maintaining a very good agreement 
with the optical spectrum. The uncertainty on log g from the optical 
analysis alone is quite large because of the modest signal-to-noise 
ratio of the X-shooter spectra. Using the FUV spectra in this case 
does help better constrain the value of the surface gravity. 

We also note that, including error bars, WLM-All is over- 
luminous (by « 0.2 - 0.3dex), compared to luminosity class I stars 
with similar spectral type in the SMC or LMC (see e.g. |Massey| 
et al.|2009^. This trend is also present, although even stronger, in 
Tramper et al.| ( |20lT||2014[ >, and is consistent with the “a” subscript 
on its classification. 

For the photospheric parameters listed in Table|^ the iron lines 
in the FUV spectra of IC 1613-Bll and WLM-All are well re¬ 
produced with models having Fe/Fe0 = 0.14 and = 17 km s“^ 
and 13 km s“^, respectively. However, minor changes to the micro¬ 
turbulence in the models (a decrease by 3 km s“^) to increase Fe/H 
to SMC-like value provides as good a fit as well. Figure|^illustrates 
the impact of such a minor variation of on the profiles of Fe IV 
and Fe V photospheric lines. Although the SNR of the COS spec¬ 
tra do not allow to mle out the lower metallicity (Fe/Fe0 = 0.14) 
option, we found that a SMC-like silicon abundance is more sup¬ 
ported by the lines of Si III (see Sect.|^ present in these late-type O 
supergiants. The synthetic spectra of IC 1613-Bll and WLM-All 
presented in the Appendix were obtained with models having an 
iron content Fe/Fe0 = 0.2, which we adopted as the metallicity of 
IC1613-B11 and WLM-All. 

As expected from their evolved status, the surface abundances 
of IC 1613-Bll and WLM-All present signs of CN processing. 
Carbon is depleted while we measured nitrogen enhancement with 
respect to the adopted baseline values for the metallicities of IC 
1613 and WLM. On the other hand, the scarcity of sensitive diag¬ 
nostics for oxygen abundance in the FUV makes its measure ques¬ 
tionable for these two stars. In the optical, several lines of O II are 
present, in principle, but are very weak and are not useful given 
the moderate signal-to-noise ratios of the X-shooter spectra. We 
checked that a simple scaling to 0.2 the solar oxygen content gives 
the best agreement with the observed spectra however. With this 
scaling, the sum of CNO scaled to 0.2 (as for the SMC) its so¬ 
lar value is also roughly conserved. Finally, we find that helium is 
slightly enriched in WLM-Al 1, while no change to He/H is needed 
to reproduce He line strengths in IC1613-B1L 

Here again, the most salient discrepancies with the study by 
[Tramper et al.|pOT4) are the wind parameters. The high mass-loss 
rates they derived are firmly ruled out, producing much stronger 
P Cygni profiles than observed (cf. Fig. |^. In our modeling, the 
(clumped) mass-loss rates are 1.5x10“® M0yr“^and 1.1x10“® 
Mq yr“^ for IC 1613-Bl 1 and WLM-Al 1, respectively. 

With these models, we could obtain a reasonable agreement 
with the observed P Cygni profile of C IV resonance doublet, the 
Si IV resonance doublet as well as N V AA1238, 1242, although 
none of our models could fit these three profiles with the same 
level of quality simultaneously. The case of IC1613-B11 is the most 



1616 1622 1628 1634 1640 1646 


Wovelength (A) 

Figure 4. In black, a portion of observed spectrum of WLM-Al 1 showing 
several lines of Fe IV around He II A1640. The sensitivity of the Fe IV lines 
tp the micro-turbulence is shown, for models having = 7 kms“^ (red), 
10 kms“^ (green) and 12 kms“^(blue). The top panel presents models 
computed with Fe/Fe 0 = 0.2 (SMC-like), the models in the center panel 
have Fe/Fe 0 = 0.14. The bottom panel shows a models with Fe/Fe 0 = 
0.2 and = 10 kms“^ (red), and a model with Fe/Fe 0 = 0.14 and = 
13 kms“l (blue). 


problematic as we could never reach an excellent fit of the C IV pro¬ 
file for this star with our models. At low velocities, the absorption 
component of the observed P Cygni profile is stronger than pre¬ 
dicted by the models, unless we increase M to such extent that in 
the synthetic profile, the absorption is too strong near Voo and we 
also produce a strong emission component that is not observed in 
C IV. Note that this behavior is typical of late-type O dwarfs ex¬ 
hibiting the "weak-wind” problem, although IC 1613-Bll is above 
the luminosity threshold usually found for these stars (log L/Lq~ 
5.2, e.g. [Martins et al.|2005[|Marcolino et al.|2009| l. Furthermore, a 
well-defined P Cygni profile is produced for the Si IV doublet, also 
not observed. 

The P Cygni profiles, although relatively weak, presented 
good sensitivity to variations of the clumping filling-factor. The fill¬ 
ing factor is moderate in both cases, that is foo =0.1, implying a 
reduction of a factor of three for M compared to an homogeneous 
case. The M we derive correspond to a regime where Ha is mostly 
insensitive to constrain mass-loss rates ( [Marcolino et al.|2009[ >. The 
modeling of the FUV and optical spectrum of supergiant stars in the 
SMC (Bouret et al., in preparation) yields mass-loss rates in good 
qualitative and quantitative consistency with those derived here for 
IC 1613-Bll and WLM-All. The mass-loss rate derived for AzV 
327, for instance, is roughly three times as high as those of IC 1613- 
Bll and WLM-All, a trend that is anticipated from the relative 
strength of their P Cygni profiles (cf. Fig.[TJ. 
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Figure 5. In black, line profiles for different transitions in the FUV range, observed for IC1613-B11 (left) and WLM-Al l(right). Overplotted in red, P Cygni 
profiles produced by models using the par ameters from|Tramper et ar]j2014). 


The exponent of the velocity law was constrained from the 
weak emission component of the P-Cygni profile. Its value is con¬ 
sistent with theoretical predictions JMuijrese^^ 201^k but is 
lower than found by |Mokiem et al.| ( |2006| > or |Massey et al.H2009^ 
for late-type supergiants. Furthermore, the wind terminal veloci¬ 
ties are consistent with Voo measured by |Massey et al.| ( |2009^ from 
UV data for late-type supergiants in the Magellanic Clouds. How¬ 
ever, we stress that the P Cygni profiles in B11 and All are weak 
enough that the wind terminal velocities are rather difficult to mea¬ 
sure. Most particularly, several lines of Fe IV overlap the blue edge 
of the C IV resonance doublet P Cygni profiles. These lines were 
not an issue for IC 1613-A13 but they are stronger at the lower Tgjy 
ofBll and All. With A V si 1200 km s“^, they fall right in the typ¬ 
ical range for Voo indicated for B11 and Al 1 by the discrete absorp¬ 
tion components, following |Prinja| ( |T998| l. Although the CMFGEN 
models do account for these Fe IV lines, their sensitivity to Tgjy, the 
adopted Fe/H, the microturbulence and even log g could introduce 
an additional uncertainty in the determination of Voo , with expected 
repercussions on the predicted mass-loss rates (see Sect.|10^. 


6 METALLICITY OF THE MASSIVE STARS OF IC 1613 
AND WLM 

The striking conclusion by |Tramper et al.|p011| > of a discrepancy 
between observed and predicted mass-loss rates relied for a large 
part on the adopted metallicities for IC 1613, WLM and NGC 3109. 
The (low) metal content for these galaxies has been measured by 
various techniques, which all point towards a SMC-type metallic- 
ity or even below. In their more recent work,[TYamper et S^j2014^ 


acknowledged that the metallicities they considered in their initial 
study were likely underestimated. In the present study, we have a 
direct access to the global metallicity of the stars through thou¬ 
sands of photospheric lines of iron ions (Fe III to Fe VI) present in 
the FUV spectral range. Oxygen abundance diagnostic is also avail¬ 
able (cf. Sect. |4.1^ . Before comparing empirical mass-loss rates as 
derived from the modeling of FUV spectra to theoretical ones, we 
first need to discuss how the metal content we derived compares to 
previous determinations. 


6.1 Oxygen 

As discussed in Sect. and |5.2| the oxygen abundance is diffi¬ 
cult to constrain accurately in the three stars. The best we can say 
at this point is that oxygen abundances such that the total CNO 
is conserved do provide good fits to the scarce diagnostics avail¬ 
able. In principle, the evolved status of IC1613-B11 and WLM- 
Al 1 should translate into oxygen depletion, compared to the initial 
CNO mixture. Although we do not need such a depletion to repro¬ 
duce the O IV lines present in the FUV spectra of both stars, we 
cannot rule out that oxygen abundances are actually lower. Only 
for IC 1613-A13 does such a depleted oxygen content manifest it¬ 
self (maybe a consequence of the earlier spectral type, hence higher 
mass, implying faster evolution), but is however moderate. In sum¬ 
mary, O/H in the three stars ranges between 1/7 - 1/5 the solar oxy¬ 
gen content, but we emphasize again that these estimates are poorly 
constrained and may be revised provided stronger diagnotiscs be¬ 
come available. How relevant these abundances are when dealing 
with comparisons between observed and predicted mass-loss rates 
of massive stars is however questionable. 
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6.2 Stellar Iron Content 

Although the low oxygen (and more generally a-elements) content 
of IC 1613 and WLM is confirmed from FUV spectroscopy, di¬ 
rect measures of the iron abundance are more relevant to study the 
metallicity dependance of mass-loss. It is expected that the mass- 
loss rates of massive stars are determined by the the radiative ac¬ 
celeration below the sonic/critical point in the wind, which is pre¬ 
dicted to be driven by the iron group elements for O-type stars in the 
metallicity regime of our sample ( |Vink et al.|2001^ . Therefore, such 
measures of the iron abundance as those we achieve in this work 
are mandatory to compare empirical mass-loss rates to theoretical 
predictions, and in any case are more relevant than metallicities 
based on a-elements abundances. Iron abundance determination of 
course relies on the strength of iron lines which depends on the ac¬ 
tual iron ionization. The later depends on the effective temperature 
and gravity (e.g. [Heap et ^|2006[ l. For a given temperature, the 
ionization is higher at lower gravities as expected from the Saha 
equation and the numerous iron lines grow stronger. Iron lines may 
also look stronger because of higher microturbulence as shown 
in |Flillier et al.| ( [200^ and |Bouret et al.| (2003) . This degeneracy 
of the ionization balance of iron with the photospheric parameters 
is such that deriving the iron abundance from fits to the iron forest 
requires Tgjy, log g and to be determined from other diagnostics, 
which we have at our disposal through the analysis of the optical 
spectra for instance. From the joint analysis of the FUV and opti¬ 
cal data, we found that for the three stars presented here, [Fe/H] = 
-0.66, i.e. their iron content is SMC-like. 

In IC 1613, stellar iron abundances of [Fe/H] = -0.67 have 
been inferred from spectroscopy but for late-type, evolved stars 
(e.g. M type supergiants [Tautvaisiene et al.||2007'] ). Alternatively, 
photometric studies were performed to measure the iron to hydro¬ 
gen ratio of stars on the red giant branch ( [Cole et al.|1999[[Skillman| 
|et al.|[2Q0^ [Minniti & Zijlstr^|1997t |McConnachie et al.||2005| >T^ 

These measures, however, determine the iron content of an older 
stellar population, while the objects we are dealing with here be¬ 
long to a young stellar population. To our knowledge, no other 
measures exist of the present-day iron abundance as probed here 
through high spectral-resolution FUV spectroscopy and quantita¬ 
tive analysis of the forest of iron lines. The iron content we derive 
for the two stars in IC 1613 is also compatible with the study of the 
star formation history of IC 1613 by |Skillman et al.| ( |2003^ ; these 
authors indeed found that the metallicity increased from [Fe/H] = - 

1.3 to -0.7 over the age of the galaxy. It is also noteworthy that using 
oxygen + silicon abundances as a proxy for the global abundance 
of a-elements, [a/Fe] is solar in our case (or slightly lower for IC 
1613-A13). Given the uncertainties in both the abundances for a- 
elements as well as iron, this is however in qualitative agreement 
with the conclusions of |Tautvaisiene et^ ( |2007) for IC 1613. Fi¬ 
nally, we note that from a qualitative comparison of FUV spectra of 
OB stars in IC 1613 with OB stars in the SMC, [Garcia et al.]j2014^ 
also argued that the iron contents are similar in both galaxies. 

On the other hand, |Urbaneja et ^ ( [2008^ quote a (present-day) 
metallicity of log (Z/Z 0 )=-O .8 (±0.2) for WLM-Al 1, although it is 
not derived from the direct analysis of iron-group lines. Although 
this value for Z/Z© agrees with our estimate of the metallicity 
within the error bars, the average metallicity of the A and B su¬ 
pergiants analyzed by jUrbaneja et al.] j2008| > is log (Z/Z0)=-O.87 
(±0.06). This lower value is not favored by the spectroscopic anal¬ 
ysis in our case, although better signal-to-noise ratio of the spec¬ 
trum of WLM-Al 1 would be required to firmly discard this so¬ 
lution. We note that the metallicity of A-supergiants analyzed in 


WLM by jVenn et al.jpOO^ is however even higher than ours (with 
[Fe/H] = -0.38 ±0.20). The ranges in [Fe/H] are large enough that 
despite measure uncertainties, they hint to a real scatter in the star- 
to-star metal content in WLM. This comes as a surprise as WLM- 
Al 1 and the stars analyzed in jVenn et al.| ( |200^ are located within 
the same region of the WLM galaxy. 

Overall, our joint analysis of the FUV and optical spectra re¬ 
veals that the iron content of the three stars is similar to the SMC. In 
this sense, the theoretical mass-loss rates are expected to be higher 
than previously considered, although by a fairly limited factor of « 
1.3. This is far from being enough to account for the strong discrep¬ 
ancy discovered by jTramper et al.| ( |201 l| l and confirmed in jTramper] 
jet al.| ( |20T4l l, especially for the late supergiants. The wind terminal 
velocities are also impacted by the metallicity, although to lesser 
extent, which does not change this conclusion. 


7 BINARY STATUS 

Multiplicity can impact the mass-loss determination in several way 
(see the discussion in jTramper et al.|2011) . The issue of multiplic¬ 
ity among massive stars has become very important in the past few 
years, with results for galactic massive stars pointing to fractions 
of stars with companions ranging from 75% ^Mason et al.||200^ 
for clusters and associations in this specific case) and up to a strik¬ 
ing 91% ( [Sana et al.|2014[ l. At lower metallicity, results obtained 
within the VLT-FLAMES Tarantula Survey ( [Sana et al.|2012| and 
refs therein), point to a fraction of more than 50% binaries within 
the sample, once corrected for observational biases. 

Concerning IC 1613-A13, the SED fitting is compatible with 
the star being single; the fluxes from the COS spectrum and the 
photometry (cf. Table are well accounted for by the SED pro¬ 
duced by a single-star model as computed with CMFGEN. The 
luminosity we derive is slightly higher (within 0.2dex) than those 
we obtained for (three) stars with close spectral type and effec¬ 
tive temperature in the SMC ( jBouret et al.||201^ . However, one 
of these stars, namely MPG 324 that we used for morphological 
comparison in Sect. (see also Table and Fig. is probably 
a single-lined binary as indicated from small radial-velocity vari¬ 
ations jEvans et al.||2006) , even though its theoretical SED from 
CMFGEN is compatible with the observed one. The case of MPG 
324 suggests that we cannot firmly rule out the presence of a com¬ 
panion for IC 1613-A13 from the SED fitting approach only. Nev¬ 
ertheless, we can safely conclude that since the putative companion 
is of low mass enough that it does not show up in the observational 
SED, it is unlikely to impact the wind signatures we observe in the 
EUV spectrum. 

The effective temperature, surface gravity and luminosity of 
IC 1613-Bl 1 are compatible with those derived for stars with same 
spectral type/luminosity class in the MCs (e.g. [Massey et al.|2009| >, 
for which no signs of multiplicity are known. Much like for IC 
1613-A13, we do not see an influence from a potential companion 
on the SED of IC 1613-B11; here again, if the star is in a binary sys¬ 
tem, it must be with a companions of late type, whose luminosities 
must be low, hence, expectedly, the mass-loss rate. More signifi¬ 
cantly, the photospheric parameters of IC 1613-Bl 1 are consistent 
within the error bars, with the calibration in [Martins et al.| j2005[ > 
for an 09.5 supergiant. Although low metallicity stars are known 
to be on average hotter and more luminous than their galactic coun¬ 
terparts for a given spectral, it is well-known that theses differences 
decrease at later spectral type and practically vanishes by BO (see 
e.g. Table 9 in|Massey et al.|2005). The spectroscopic mass is also 
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consistent with this calibration, which further supports the single 
star interpretation. 

The case of WLM-All is somewhat different. The SED fit¬ 
ting indicates that this star is more luminous than a normal 09.7 
supergiant (also denoted by the ”a” subscript of its spectral classifi¬ 
cation). The photospheric radius we derive and the surface gravity 
we measure are such that the star is then roughly twice as mas¬ 
sive as other late-type supergiants of similar spectral properties, in¬ 
cluding IC 1613-Bll or stars in |Massey et al.| ( |2009| l. Such a high 
spectroscopic mass (^53 Mq) is actually more characteristic of 
O-supergiants with early spectral types (e.g. |Massey et al.||2009| 
[Bouret et al.|2012^ . This might indicate that WLM-Al 1 is actually 
in a binary system or that there is a unresolved star in the fore/back¬ 
ground, contributing to the observed fluxes. We will come back on 
the influence of changing the stellar luminosity on mass-loss pre¬ 
diction, in a following section. 


8 FUV MASS-LOSS RATES AND UNCERTAINTIES 

[Tramper etal.| ( |20lT||2014| l relied on the fitting of the Ha (and He II 
A4686 line for the hotter objects) to measure the mass-loss of their 
targets. Known advantages are an easy and robust abundance deter¬ 
mination and ionisation balance calculation. The sensitivity of Ha 
is limited to stronger winds however, showing only moderate line 
core filling for mass loss rates between 10 “® and 10 “^ Mq yr“^ 
(see Fig. 5. in jMarcolino et al.|2009) . Because of their low metallic- 
ity, our targets should have mass loss rates in this parameter range 
( jVink et al.|2001^ . Detailed investigations should therefore rely on 
diagnostics showing sensitivity in this regime, namely FUV. Ob¬ 
viously, mass loss determinations in this spectral range also suffer 
from specific uncertainties related to the assumptions adopted in 
the modeling. 


8.1 Abundances, Ionization and Wind Acceleration 

The observed COS spectra present three important lines of carbon 
ions, namely C IV A1169, C III A1176 and C III A1247, that we 
used to measure the photospheric abundance of carbon. Within the 
uncertainties on C/H, the mass-loss rates of the targets could vary 
by at most a factor of two (but see below for a comment on the influ¬ 
ence of M on the C III Al 176 line in the late-type supergiants). We 
emphasize that the predicted (wind) profiles of the Si IV resonance 
doublet are reasonably reproduced in our models of the late super¬ 
giants (cf. Fig. |A3| and Fig |A5[ ). The Si element does not suffer, a 
priori, from the uncertainties related to carbon, as it is not involved 
in nuclear burning until very much later phases of the evolution 
of massive stars. Changing M by a factor of two (i.e. within the 
boundaries of the relative uncertainty due to the carbon abundance) 
would change the profiles of Si IV resonance doublet enough that 
we can rule out these solutions. 

The other obvious uncertainty is related to the ionization frac¬ 
tion of CNO ions. This might be a potential flaw of the FUV- 
based mass-loss determination, compared to hydrogen-helium 
based methods for photosphere -|- wind properties determination. 
This is because of uncertainties in the predicted flux levels at 
200 sC As: 400 A, or more exactly differences in the flux levels 
predicted by state-of-the-art NLTE, line-blanketed models for mas¬ 
sive stars ( jPuls et al.|2005| >. These differences translate in differ¬ 
ences in the predicted populations of important ions such as C IV 
or N V for instance, hence on the estimate of the mass-loss rate. 


[Puls et ar| ( |2005| l emphasized that the uncertainties on the flux lev¬ 
els mostly impact the ionization fractions in the outer part of the 
winds, where they translates into differences in population of typ¬ 
ically ± 0.15 dex, while they agree almost perfectly everywhere 
else. In our case, we are able to reproduce line ratios of differ¬ 
ent ionization stages of carbon, nitrogen, and oxygen which should 
guarantee that the flux slope in the critical range 200 A ^ 400 
A is correct. In any case, the ionization fraction of carbon, nitrogen 
and silicium ions should be accurate within a factor of a two, which 
does not change our conclusion on the actual mass-loss rates of our 
targets. 

We can obtain a good grip on the wind velocity law, accel¬ 
eration parameter and terminal velocity, from fitting of the slope 
of the emission component of the P Cygni profiles. It is in any 
case more robust than what can be done from optical lines alone 
when the line profiles are in absorption, in which case these pa¬ 
rameters are not accessible thus increasing the uncertainty in the 
derived mass-loss rate, by a factor from 2 to 3 ( jPuls et al.|1996^ . In 
our case, the uncertainty in /3 and Uoo reflects in a moderate cor¬ 
rection on M only, which again does not change the conclusion of 
the present paper. Note that uncertainties on Voa have an influence 
when dealing with the modified wind-momentum-luminosity rela¬ 
tion (see e.g. [Kudritzki & Puls|2000|. The values of I4o we mea¬ 
sure are for instance « 25% different from those listed in|Tramper| 
|erari([20T4ll, higher for IC1613-A13, but smaller for IC1613-B11 
and WLM-All. 


8.2 Clumping 


Wind clumping and the formalism used to describe it, is another 
important issue for the determination of mass-loss rates. It modifies 
the predicted line profiles for a given M and this parameter cannot 
be considered separately of clumping parameters anymore. 

Clumping as implemented in CMFGEN is usually referred 
to as micro-clumping jHillier & Miller|1999| and refs therein). Its 
influence on M determinations in O stars, both from FUV and op¬ 
tical diagnostics has been tested and discussed extensively for stars 
in the Galaxy and the SMC (e.g. [Hillier et al.|2003| [Bouret et al.j 
[2005[[20T2][2013[ |. The volume filling factors (/oo) associated with 
the clumped models for our three targets (/oo = 0.03 for IC16-A13, 
and /oo = 0.1 for IC1613-B11 and WLM-All) are in line with 
those found in these studies. They imply that the mass-loss rates 
are scaled down by factors 6 and 3 (respectively), compared to the 
homogeneous values that would be derived from fitting Ha. Note 
that with the adopted parameterization of clumping inCMFGEN 
(see Sect. |^, the volume filling factor in the formation region for 
the Ha line is much higher than the value /oo. 

However, the filling factor approach is admittedly too simple. 
It considers optically thin clumps, an assumption valid for the con¬ 
tinuum but not so much for the lines. Furthermore, it does not ac¬ 
count for the interclump medium which is important for producing 
the P Cygni profiles of super ionized species jZsargo et al.|2008) > 
or, even more fundamen tally, because void interclump medium is 
physically unlikely (e.g. 


Owocki et 


ar]|l988 


Surlan et 


'ZII2OI3 


[Sundqvist et al.|2014[ >. It also ignores porosity, both in spatial and 
in velocity space (e.g.|Oskinova et al.|2007[ [Sundqvist et al.|201 1[ (. 

A detailed study by [Sundqvist et al. (2014 1 concluded that velocity- 
porosity, a.k.a vorosity, impacts both the line profiles (diagnostics) 
and mass-lo ss properties (dynam ics), in contradiction with an ear¬ 
lier result by Surlan et al. (2013 1 and based on a different approach, 
showing that spatial porosity was most important. 

Notwithstanding the very different levels of physical sophisti- 


(gl 2013 RAS, MNRAS 000,pT[^ 















































13 


cation to describe the nature of wind clumping and its effets (e.g. 
the role of porosity/vorosity), studies have found that the correc¬ 
tions to mass-loss with respect to theoretical predictions (the widely 
used recipe by | Vink et al. 120011) should be less than a factor of 1.5 


2.5 typically ( Bouret et al.|2012| 


Surlan et al.l2013 


Sundqvist 


|et al.|20l4l l. This difference is actually very reasonable given the 
uncertainties related to both approaches, namely NLTE radiative 
transfer codes used to describe clumped stellar winds on one hand, 
and Monte-Carlo calculations of the radiative force used to predict 
mass-loss rates on the other hand. 

Although these results were obtained for the same early-type 
Galactic supergiants with fairly dense winds, it is tempting to ar¬ 
gue that the mass-loss rate we obtained, at least for IC1613-A13 
which has a relatively dense wind, is realistic within the same mar¬ 
gin, enough to provide a robust comparison basis against theoretical 
predictions. 


8.3 Line force, momentum equation and implications 

Whatever the physical nature of clumping, the presence of inhomo¬ 
geneities in the wind makes it very difficulty to determine the sys¬ 
tematic errors associated with our mass-loss rates estimates. This 
is dictated by the lack of diagnostics, a problem exacerbated in 
stars that have weaker winds. Typically in the FUV we only have 
a few lines suitable for mass-loss determinations, and these line 
often have other sensitivities. For example in A13, the only wind 
diagnostics are the C IV and N V doublets, N IV A1718, and the ab¬ 
sence of O V A1371, while in All we loose O V as a diagnostic 
but gain the Si IV doublet. For an inhomogeneous wind we have 
multiple parameters that describe the the clumping, porosity, and 
porosity in velocity space that need to be determined. 

In our modeling we have adopted the volume filling factor ap¬ 
proach which has only two free-parameters, and one of which (the 
location at which clumping starts) is not well constrained. How¬ 
ever, we have one additional consistency check that can be utilized. 
While CMFGFN does not compute the mass-loss rate self consis¬ 
tently, it does compute the radiative force, hence we can check how 
well the momentum equation is satisfied. 

For A13 we find that the momentum equation is not satis¬ 
fied - we deposit too much momentum in most of the wind, and 
this inconsistency cannot be changed by altering^the velocity law 
while remaining constant with the observation^ Three possible 
solutions to the inconsistency are: 


(i) An incorrect ionization structure. 

(ii) Additional hot gas that we are not detecting. In the standard 
wind-driven wind model, it is assumed that hot shocked gas, aris¬ 
ing from the radiation driven wind instability, rapidly cools. While 
this is likely to be valid for the wind of a star like C, Pup (see, e.g., 
discussion by |Hillier et al.|[T993) , the assumption will eventually 
become invalid as the density declines jDrew et al.||1994[ |Mar-| 
|tins et al.|2005[[Lucy|2012) . Recently |Huenemoerder et ~ j20 1 2) 
showed that the X-ray emission from ^ Columbae originated in the 


^ We ignore issues in accelerating the wind around the sonic point. In gen¬ 
eral we have too little line force at, and just above, the sonic point to accel¬ 
erate the flow. Much of this discrepancy arises from issues with the choice 
of the microturbulent velocity: driving a stationary wind through the sonic 
point requires the microturbulent velocity to be much less than the sound 
speed. Although this is typical for O supergiants, it is not the case here - 
the microturbulent velocity is of order 0.5 times the sound speed. 


wind, and required a hot gas component with a mass-loss rate in 
excess of that inferred from analysis of UV lines. 

(iii) We are underestimating the mass-loss rate due to the effects 
of porosity and/or vorosity. 

We dismiss the first one since in CMFGFN the predicted and 
observed ionization are consistent (see comments above). For AI3 
we also argue that in its relatively dense wind, the hot gas com¬ 
ponent is probably not important. However, consistent with other 
findings, porosity and vorosity could be important. If we adopt the 
same /oo and M=5.5x 10“^ Mq yr“^ we find much better agree¬ 
ment with the radiation force and that required to derive the wind. 
However, absorption features in the spectrum are now too strong. 
For /oo = O.I and M = 4.0 x 10“^ Mq yr“^ the situation is some¬ 
what more intermediate - we somewhat deposit too much momen¬ 
tum between 100 and 1600 kms“^, and somewhat too little in the 
outermost regions (V > 2000km s“^). 

Thus for A13 we conclude that our derived mass-loss rate is 
underestimated, at most, by a factor of two. This is consistent with 
our finding in (^Pup where we found M= 2.0 x 10“® Mq yr“^ 
for fao = 0.05 using multiple diagnostics. However we could only 
fit the P V resonance doublet with an abundance a factor of two 
less than solar. If we assume the solar abundance is applicable, 
this would imply that vorosity/porosity effects influence mass-loss 
rates derived from the P V doublet by about a factor of two. 

For All and Bll, with much weaker winds, the situation is 
much more problematical. Let us discuss the case of All as an 
example. The mass-loss rate of 1.1 x 10“® Mq yr“^ is a factor of 
three less than the limit for a single optically thick line located at 
the emission peak of the SED {Ljc?). This inconsistency is borne 
out when one examines the momentum equation - the radiation 
force that we derive is typically a factor of 5 to 10 times large than 
that required to drive the flow. This discrepancy indicates that some 
fundamental assumption in our modeling is incorrect. 

For All we believe that the problem is twofold. In its weak 
wind it is possible that the gas does not cool after it is shocked. 
Thus a significant fraction of the wind is in hot gas - a result con¬ 
sistent with the use of X-ray filling factors near unity in order to get 
a logLx/Lboi « — 7. The second is that in these weaker winds, the 
influence of porosity/vorosity is more important. An alternative is 
that wind dynamics is influenced by the presence of a strong mag¬ 
netic fields but since the issue of weak winds in persistent, even in 
stars where a magnetic field is not detected (e.g., 10 Lac, |David-| 
lUraz et al.|20l4l ), this cannot be the answer. If we adopt f^a =0.1 
and M = 3 X 10~^ Mq yr“^, there is reasonable consistency be¬ 
tween the radiative force and that required to drive the flow. Of 
course this mass-loss rate produces a spectrum very inconsistent 
with that observed. 

Independent of the cause of the discrepancy, we might expect 
our assumptions to be more valid deeper in the wind. Thus we de¬ 
cided to consider a model in which we varied the mass-loss rate as 
a function of the radius r, thus mimicking a wind in which, for ex¬ 
ample, a lot of the gas becomes hot and does not contribute to the 
observed UV spectrum. In the test case we assumed M= 7.7 x 10“® 
Mq yr~^ in the inner wind, and 1.1 x 10“® Mq yr~^ in the outer 
wind, with a smooth transition between them. With this model two 
distinct features were revealed: 

(i) The higher mass-loss-rate model desaturates the C III triplet 
line near 1176 A, and produces a weak blue wing. The profile is in 
better agreement with observations, than our best fit model. 

(ii) The higher mass-loss-rate model produces a slight blueshift. 
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Figure 6. Variation of the mass-loss rate, plotted here as a function of the 
velocity in the wind. The model was used to compute the spectrum pre¬ 
sented in Fig. ^(blue line). 




Figure 7. Profiles of C III A1176 (top) and Si IV resonance doublet (bot¬ 
tom). The model with M varying with radius (see text and Fig.|^ is plotted 
in blue. Our best fit model is shown in red. 


and a slight broadening, of the absorption associated with the Si IV 
resonance doublet. This blueshift and broadening is also seen in the 
C in triplet. 

We thus argue the our mass-loss rate for All derived from 
the conventional lines is probably too low, by some factor between 
3 and 10. Again, similar conclusions hold for Bll. In these weak 
wind stars, the influence of porosity/vorosity effects, and the possi¬ 
ble presence of hot gas, can potentially severely bias the estimated 
mass-loss rate. 


9 WIND PROPERTIES AT LOW METALLICITY 


From the modeling of the optical spectra of early-type stars, in 
three metal-deficient galaxies, more precisely galaxies with O/H 
less than SMC, |Tramper et al.| ( |2011| |2014) strikingly concluded 
that the mass-loss rates are consistent with those that would be pre¬ 
dicted by the radiatively-driven wind theory for a higher, LMC- 
type, metal abundance. This lead them to argue that this theory 
would break down at sub-SMC metallicity. If confirmed, this would 
be a major problem for our understanding of the physics of stellar 
winds, and for stellar evolution in the upper HRD as well. How¬ 
ever, several facts suggest that the conclusion by [Tramper et al.| 
(|201 1[[20T4) may be premature. 


• We showed in Sect.|^that the iron abundance in IC 1613 and 
WLM is higher than the metallicity adopted from O/H ( [Tramper] 
|et al.|20rT| and refs therein). In the three stars, the metal content 
is not, if at all, significantly smaller than in the SMC, f or which 


theoretical M were found to match the empirical ones (Mokiem 


|et al.||2007] ) . The higher metal content shall impact the theoreti¬ 
cal mass-loss rates as well as the wind terminal velocities, hence 
the modified wind-momentum-luminosity relation (WLR) for the 
star0 

• No important change in wind driving (e.g. going from a driv¬ 
ing through lines of iron group elements to CNO) is expected until 
the metallicity is I/SOZq and below ( [Vink et al.|2001 ^. 

• [Herrero et al.| \2Q12) found that the wind parameters of an 
06.5 Illf star in IC 1613 are such that the star follows the WLR for 
its metallicity. In other words, if the metallicity would be sub-SMC 
(but see the recent conclusions by |Garcia et al.|20T4l besides ours), 
there is at least one case where the standard theory still applies. 

• The sample analyzed by |Tramper et al.| j2014| > contains 10 
stars, a majority of them having wind strengths that are compatible 
with theory within the errors, assuming a SMC-like metal content. 
Only three stars (all of them in NGC 3109) are real outliers in the 
WLR diagram. 

• The intrinsic dispersion in the wind parameters within a given 
spectral/luminosity class typically ranges within a factor of two. 
This is of the same amplitude as the difference that is expected for 
two stars with the same parameters but a factor of 2.5 in the metal 
content (i.e. going from LMC to SMC-type metallicity). 


We showed in Sect. |5.1| and |5.2| that the FUV-spectra of IC 
I6I3-AI3, IC 16I3-BII and WLM-AII indicate the presence 
of significant wind clumping. All other things being equal, the 
mass-loss rates quoted in |Tramper et al.|p014| l should therefore 
be reduced accordingly. The amplitude of the reduction depends 
on the nature of clumping (cf. Sect. |8.2^ ; the typical reduction of 
M found for clumped versus homogeneous models is found to be 
a factor 3 to 6 . Clumped mass-loss rates are reduced by less than 
a factor of three cor npared to theoretical pre dictions ^Bouret et al.| 


2012 


Sudan et al. 


2013 


Sundqvist et al. 


2014 I. Note however 


that this last conclusion holds as long as the impact of clumping 
on the wind driving is not considered (see e.g. |Mu ij res et al.| 
|201 1| for comments). Recently, [Sundqvist et al.H2014 i found that 
when applying a vorosity correction to the line force (based on 
the standard theory of line-driven winds by [Castor et al.||1975) >, 
a reduction of the line force that drives the wind is obtained. 


® This relation links the modified wind momentum, defined as Dmom = 
M Voo v/ R/Rq, to the stellar luminosity. The WLR is often used to study 
mass-loss in various galaxies and environments (e.g. [Mokiem et al.|2007[ 
as Dmom is almost independent of mass. 
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Therefore, a reduction of theoretical mass-loss rates is expected 
from these calculations if vorosity is present near the critical point, 
where the mass-loss rate is set. Such a scenario is supported by 
several observational (e.g. |Puls et ^ 2006 Najarroetay|20]T 


Bouretetay201^ and theoretical ( Cantiello et al.|200^ Sundqvist| 
& Owocki [2013^ results suggesting that clumping indeed starts 


close to the photosphere. 


The mass-loss rates we derived using models with clumping 
are lower than those obtained from the modeling of the optical 
spectra only, especially for the late-type supergiants. We present 
in Fig.[^the synthetic profiles of wind sensitive lines in the optical, 
namely the He II A4686 and Ha line, as predicted with our values 
for M and those obtained with CMFGEN for the parameters listed 
in ( [Tramper et al. [20141 their tables 4 and 5). 

The Hq line in IC1613-A13 are very similar in shape and 
strength, but a striking difference is observed in the predicted pro¬ 
files for He II A4686. Even more striking is the fact this line is pre¬ 
dicted in absorption by FASTWIND jTramper et al.|2014l cf. fig¬ 
ure 4 in), while the profile predicted by CMFGEN for the same pa¬ 
rameters (and no clumping) shows clear electron scattering wings 
in emission. It is known however that CMFGEN often produces 
more emission in Ha and He II A4686 than FASTWIND, for in¬ 
stance in the temperature and surface gravity regime corresponding 
to IC16I3-AI3 (e.g. |Rivero Gonzalez et al.|2QI2^ . In this regime, 
the formation of He II A4686 is very sensitive to radiative trans¬ 
fer effects, for instance the detailed treatment and amount of line¬ 
blanketing (model in green in Fig. has much less elements than 
the best fit models), and the line profile responds accordingly. 

Concerning the late-type supergiants IC1613-B11 and WLM- 
A11, the profiles predicted by CMFGEN for He II A4686 and Ha 
for the two sets of parameters bracket the observed profiles. Since, 
admittedly, CMFGEN predicts stronger emission than EAST- 
WIND in these lines, these profiles in Fig. also likely bracket 
those which lead [Tramper et ar] ( |2011[|2014[ rto claim that mass- 
loss rates are higher than expected in IC1613 and WLM. We argue 
that, from such profiles only, it is quite uncertain what M should 
be favored. 

Despite the uncertainties on the real absolute values of mass- 
loss rates introduced by the feedback from clumping on line driv¬ 
ing, the considerations spelled out hereabove suggest that the em¬ 
pirical mass-loss rates derived by |Tramper et al.| ( |2014^ are likely 
overestimated, compared to those expected theoretical predictions 
and those derived by clumped models applied to the FUV -l- optical 
spectral range. 


10 OBSERVED VERSUS PREDICTED MASS-LOSS 
RATES 

10.1 Semi-empirical mass-loss rates 

The mass-loss rates calculated for the stellar parameters using the 
recipe by |Vink et al.| ( |200T] l are presented in TableFor the three 
stars, M measured from FUV spectroscopy are lower than these 
estimates. 

In the case of A13, the difference is slightly more than a factor 
of three, which is reasonable given the uncertainties on the photo- 
spheric and wind parameters derived from spectroscopy. Note that 
this difference would be sightly reduced if Z = 0.14Z0 applies in¬ 
stead to the stars (by a factor (0.2/0.14)°'®®“°'®® = 1.28-1.34), al¬ 
though the spectroscopic analysis does not support this lower value. 


The mass-loss rate and clumping parameters of A13 are consistent 
with recent results obtained for SMC dwarfs ( [Bouret et al.|2013l >. 
The same behavior between measured and theoretical mass-loss 
rates is obtained for early type galactic dwarfs jBouret et al.|2005| 
[Martins et al.|2012^ and supergiants ( [Bouret et ^.[2012 1 . All in all, 
this indicates that the wind properties of O stars of early spectral 
type are well accounted for by the standard theory of radiatively 
driven winds, regardless of the luminosity class and stellar metal- 
licity. As for the wind terminal velocity, we note that our value V/o 
= 2180 kms“® is much higher than 14o ~ 1570 kms“®, inferred 
from the usual scaling Voo/Vesc = 2.6 ([Kudritzki & Puls[2000^ , then 
to the star’s metallicity V/o oc Z°'®® ( Leitherer et al.|19f^ l. It is im¬ 
portant to recall that in the recipe by Vink et al. (2001 1 , the wind 
velocity structure and terminal velocity are adopted as input to pre¬ 
dict the mass-loss rate. Said differently, the adopted value for Voo 
impacts the mass-loss prediction. A difference in I4o such as the 
one we find for A13 implies a difference of more than 0.28 dex in 
the theoretical mass-loss rates, all other parameter being equal. 


The case for B11 and Al 1 is different. The terminal velocities 
we measure are in good agreement (within 50 km s“®) with the val¬ 
ues predicted using the theoretical scaling between surface escape 
velocity and wind terminal velocity, then with metallicity (for the 
derived stellar parameters and iron abundance). On the other hand, 
the mass-loss rates are much lower than predicted by the recipe 
from [Vink et SI] ( [200 1^ for these stellar parameters and the iron 
abundance. The discrepancy reaches up to a factor of twenty and 
thirty five for Bll and All (resp.), assuming a SMC-type metal¬ 
licity for IC 1613 and WLM. In order to reconcile the measured 
mass-loss rates with theoretical values, the metallicities would need 
to be Z = 0.03 Zq and O.OO2Z0 for Bll and All respectively. Such 
low metallicities are clearly ruled out by the spectroscopic analysis 
(as is Z = O.14Z0, which would only marginally reduce the dis¬ 
crepancy between observed and theoretical M anyway). Of course, 
this discrepancy is very likely an upper limit, should the true mass- 
loss rates be revised upward for All and B1 las a result of a proper 
treatment of porosity/vorosity effects (cf. Sect.[^. Note also that in 
the specific case of All, the difference could be reduced though, 
adopting a luminosity and a mass more in line with those expected 
for a normal 09.7 supergiant, to compute the theoretical mass-los s 
rate (e.g. log ^ = 5.5 and M=31M0, from 


Martins et al. 


2005 I. 


When corrected for this luminosity effect, possibly indicating a bi- 
nary/multiple status for this star, the difference with the theoretical 
Mvink is of the same amplitude as for B11. 


In any case, even though our modeling concurs with [Tramper[ 
[et al.[ ( |^ll[[2014) that the measured mass-loss rates of Bll and 
All disagree with theoretical predictions, we find the discrepancy 
to occur in the opposite sense. 


Recently, [Muijres et al.[ \20\2\ built on the work by [Muller[ 
[& Vink[ ( [2008) to compute new mass-loss rates and terminal wind 
velocities, improving the physics implemented in the Monte-Carlo 
(MC) method used by [Vink et aT| ( [2000[ [2001[ l. In this new ap¬ 
proach, a parametrized description for the line acceleration is used 
to derive new expression for the velocity field of a radiation-driven 
wind. A major conclusion of this work is that for stars more lumi¬ 
nous than log ^ ^ 5.2, the new mass-loss rates and those by 
[Vink et ^ ( [2000) agree well. On the other hand, for physical pa¬ 
rameters relevant to a typical late-type O supergiant (i.e. close to 
those of Bll and All), the new mass-loss rates are between 0.3 
- 0.4 dex lower than those computed by [Vink et ^ j2000[ l. The 
wind terminal velocities however, are notably higher than observed 
for such objects (e.g. [Howarth & Prinj^[1989[l. Notwithstanding 
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Figure 8. Left: In black, the spectral region around the He II A4686 line, observed for IC1613-A13, Bll and WLM-All(from top to bottom). Synthetic 
profiles produced by models computed with CMFGEN using the parameters from |Tramper et S!]j2014) and those derived from the joint FUV+optical 
analysis presented here are overplotted (in blue and red, respectively). Right: Same, from the region around the Ha line. The model in green presented for 
IC1613-A13 illustrates the sensitivity of the He II A4686 line to line-blanketing (Ne, Cl, Ar, Ca are not included) while the Ha line is much less affected. 


the possible reasons for these over-predicted Va^ (see|Muijres et al.| 
|2012[ for a discussion), these MC simulations essentially predict 
the total wind kinetic energy MVoo^, and any change in the wind 
terminal speed reflects in a change in the mass-loss rate, hence the 
observed decrease of M. More quantitative comparison to B11 and 
Al 1 is not possible, however, because the mass-loss rates and wind 
terminal velocities have been computed by |Muijres et al.| ( |2012^ for 
galactic metallicity only. Assuming that the theoretical scaling ob¬ 
tained by |Vink et al.H200T) still holds, the mass-loss rate should be 
scaled down by another 0.5 dex, while Voo would change by less 
than 20%, compared to the galactic values. With a total scale down 
of ~ 0.8 - 0.9 dex, M would still be a factor of three to five higher 
than measured from FUV spectroscopy. Furthermore, if the Voo for 
these late-type supergiants are much lower (as we measure them) 
than predicted (40% typically), then the theoretical M should go 
up. It is not clear that mass-loss prediction at different metallici- 
ties, based on |Muijres et al.| ( |2012^ improvement in the treatment 
of the line force parametrization, would be able to account for the 
observed mass-loss rates of B11 and All. 


10.2 Mass fluxes obtained by solving the dynamics of 
radiatively-driven winds 

|Lucy| ( (^10b^ argued that the mass-loss rates by | Vink et al.H200f) 
that have been compared to observations, may not be relevant for 
such comparisons because they are derived by imposing a global 
dynamical constraint, i.e. assuming that the energy that is needed 
to drive the wind is indeed extracted from the radiation field, rather 
than being obtained by solving the equations governing the dy¬ 
namics of radiatively-driven winds. In order to fully guarantee that 
the derived mass-loss rates are consistent with stationary transonic 
flows, |Lucy| ( |2010b^ computed mass fluxes that allow the wind 
mass flow to accelerate continuously from sub- to supersonic ve¬ 


locities. The theory of moving reversing layers (MRL) pioneered 
by |Lucy & Solom^jl970^ was used to determine the mass fluxes 
from first principles by imposing a regularity condition at the sonic 
point. An important conclusion from this work was that the so- 
called weak-wind problem mostly arises from the existence of two 
separate regimes of ionization of iron. More precisely, a shift in 
the ionization balance of iron, from Fe V to Fe IV (as Tgff de¬ 
creases), leads to weaker line acceleration below the sonic point, 
which translates into lower mass-loss rates. In the temperature do¬ 
main of late O-type stars, namely Tgjy « 30,000 K, the Fe V con¬ 
tribution to the driving is dramatically reduced and a reduction by 
up to ~ 1.4 dex is obta ined compared to mass-loss rates predicted 
by lVink et al!] ( |2000p 


More recently, LucylpOT^ extended this work to predict mass 
fluxes that depend on the stellar metallicity, without relying on scal¬ 
ing laws (e.g. |Vink et al.|200T||Mokiem et al.|2007^ but implement¬ 
ing curve-of-growth effects in the dynamic of the reversing layers 
instead. Mass fluxes computed with this method present manifest 
departure from the expected, linear, monotonic, decline with metal¬ 
licity (see Fig. 1 in |Lucy|2012^ . It is noteworthy that the strongest 
departure from a (single) power-law scaling is found for the metal¬ 
licity and (Tgff - log g) regime of the two late-type supergiants IC 
1613-Bll and WLM-Al 1. 


We used the value listed in Table 1 from |Luc^j2012^ to calcu¬ 
late the mass-loss rates corresponding to the stellar parameters we 
derive for IC1613-A13, IC1613-B11 and WLM-All respectively. 
It is fortunate that the case for T'gjy= 42,500 K and log g= 3.75 was 
also worked out by |Lucy| ( |2012[ l as it best account for the parame¬ 
ters of IC1613-A13 as derived in the present study. The results are 


® iMuijres et 31.^201^ do solve the full momentum equation, but they now 
obtain that for stars having log ^ 5.2, i.e. stars in the “weak wind’ 
regime, their hydrodynamical method fails to drive a wind. 
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summarized in Table[^ The striking conclusion from this compari¬ 
son is that, while the good match between theoretical and observed 
M for IC 1613-A13 is maintained, the discrepancy between ob¬ 
served and theoretical mass-loss rates we found for IC1613-B11 
and WLM-Al 1 when using mass-loss rates from |Vink et al.|j2001^ 
is dramatically reduced (and would be even more if mass-loss rates 
are higher by a factor of 3 to 10, as discussed in Sect [^. Quanti¬ 
tatively, the theoretical MLucy of B11 and All are more than one 
order of magnitude less than My ink - For a large part, this is a con¬ 
sequence of the aforementioned reduction of line-driving caused 
by the lower ionization of iron. The amplitude of the decrease is an 
additional balance between the decreased contribution of Fe (and 
Ni for that matter) to the driving and an enhancement of M caused 
by the lower surface gravities (see Sect. 4 and 5 in |Lucy|2010a[ for 
a discussion). The remaining differences are well within the uncer¬ 
tainties caused by the cumulative effects of error on the measured 
stellar and wind parameters (especially the impact of clumping on 
the driving) on one hand, and in input abundances (global scaling 
to solar values) and micro turbulence velocity from the TLUSTY 
models of |Lanz & Flubeny] ( |2003[ l used by |Lucy| ( |2012| l to compute 
mass fluxes. 

Interestingly, |Lucy| ( |2012[ l found no dynamically consistent 
solutions that would yield higher mass fluxes, within the purely 
radiation-driven winds framework, to provide better agreement 
with the mass-loss rates by [Tramper et al.| ( |201 l| l. An (yet un¬ 
known) additional driving mechanism would be therefore required 
to provide more radiative acceleration. With A4 derived from FUV 
spectroscopy, such a putative mechanism may not be needed any¬ 
more. 



L/Lo 


Figure 9. Modified wind-momentum luminosity relation for the sample 
stars. Black full squares mark Dmom calculated for the spectroscopic pa¬ 
rameters listed in Table|^ Blue full circles are Dmom calculated with the 
mass fluxes from |Lucy]j2012^ for the stellar parameters, while red full di¬ 
amonds are Dmom computed with the | Vink et al.|j2b01^ recipe. Dmom 
for WLM-Al 1 assuming log L / Lq = 5.5 are indicated as empty symbols 
(same shape and color coding as previously). The dashed red lines are the 
theoretical relation for LMC and SMC metallicities predicted by |Vink et al.| 
j2000) . Pink triangles are Dmom from |Tramper et aT]p014^ for the three 
stars. For clarity, no individual error bars are given but we plot a typical 
error bar on the upper left corner. 


10.3 Modified Wind-Momentum-Luminosity Relation 


Using the FUV-measured mass-loss, and those computed with 
I Vink et ^ ( |200l| l and |Lucy] ( |2012 i, we built the modified wind- 
momentum-luminosity (see e.g. |Kudritzki & Puls||2000) relation 
for the three stars (Fig. [^. We further added the modified wind- 
momenta (Dmom cf. Sect. 1^ obtained for mass-loss rates pre¬ 
dicted by both approaches for All assuming a lower luminosity 
(hence radius) for this object. Recall that this assumption is based 
on considerations on measured surface gravity, which yield abnor¬ 
mally high spectroscopic mass for this object, if the luminosity de¬ 
duced from the the SED fitting of observed fluxes is used. Mea¬ 
sured wind-momenta are clearly lower (more than one order of 
magnitude for Bll and All) than predicted using the theoretical 
coefficients ( |Vink et al.|2001| > for the WLR of the LMC and SMC 
metallicities (see caption for details). For the stellar coefficients de¬ 
rived from spectroscopy, the mass-loss rates predicted by |Vinketal.| 
( | 2001 ^ yield Dmom that are (expectedly) in better agreement with 
the theoretical relations, although some striking differences are ob¬ 
served. Note the case of IC1613-A13, where the predicted Dmom 
is more compatible with the theoretical WLR for LMC. The pre¬ 
dicted Dmom for WLM-Al 1 is below the theoretical relation but 
the agreement would be much improved assuming that WLM-Al 1 
has a log I//Z /0 lower by 0.2 dex. 

Although the small number of points makes further quantita¬ 
tive interpretation questionable, we argue that wind properties de¬ 
rived from FUV diagnostics are better accounted for by the MRL 
approach. The above comparisons strongly suggest that the infor¬ 
mation on wind quantities available in the optical spectra in [Tram- 1 
|peretal.| ( |2011[[2014| l is probably too limited and the FUV provides 
more realistic estimates, while a proper treatment of the radiatively- 
driven winds dynamics considerably improves the calculation of 


mass fluxes. The remaining discrepancies between the observed 
and theoretical mass-loss rates for A13, Bll and All reflects the 
many uncertainties related to simplifications that are made both in 
the diagnostic models and in theoretical calculation of mass-loss 
rates. The mass-loss - metallicity scaling breakdown at sub-SMC 
metallicity is not confirmed at this step. 


11 SUMMARY AND CONCLUSIONS 

We have modeled the FUV spectra of three O stars in the low metal¬ 
licity galaxies IC 1613 and WLM, observed with HST/C05. Our 
goal was to coroborate or refute the results by |Tramper et al.|j201 1| 
|2014| l suggesting that the mass-loss rates of some of these stars 
are more compatible with higher metallicity and similar to those 
of LMC stars. This result, if upheld, would have far-reaching con¬ 
sequences beyond our understanding of radiatively-driven winds, 
and would impact for instance the predicted number of collapsars 
and of SN Ib and Ic in low metallicity environments, hence in the 
high-z Universe. 

Because of numerous iron lines and lines sensitive to wind 
properties, FUV spectroscopy is the adequate tool to fully address 
and resolve this outstanding issue of the dependence of hot, massive 
star mass-loss rates with metallicity. A comparison of these new 
COS spectra with existant STIS or COS spectra of Galactic, LMC 
and SMC stars provided a direct, model-independent check of the 
mass-loss - metallicity relation. A quantitative analysis was then 
carried out using state-of-the-art NLTE unified model atmospheres 
calculated with the CMFGEN code to establish the wind properties 
of the three low-metallicity massive stars. 

• The spectroscopic analysis indicates that the metallicity of 
the three targets is I/SZq, i.e. SMC-like, although the usual value 
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adopted for IC1613 and WLM of Z = IHZq is not ruled out. Lower 
and higher values of the metal content (1/10 Zq and 1/2 Zq respec¬ 
tively) are not consistent with the observed spectra. 

• The FUV data allow secure inferences on the CNO abun¬ 
dances, of critical importance for mass-loss rates determinations. 

• The wind parameters derived from the FUV-spectroscopic 
analysis of the observed HST/COS spectra provide very satisfac¬ 
tory fits to the optical spectra as well. This further illustrates the 
limitations related to M determinations based on the optical spec¬ 
trum (Ha) only. 

• The mass-loss rates differ from those derived by|Tramper et al.| 
( |201 H[20T4l >. The difference forIC1613-A13 is moderate (less than 
a factor of 6), and it can be accounted for by clumping (and to some 
extent to the detailed treatment and amount of line-blanketing). For 
the late-type supergiants (IC1613-B11 and WLM-Al 1), the ampli¬ 
tude of the difference is much higher, reaching from ~ 1.6 to 2.5 
orders of magnitude, respectively. Here, wind clumping (absent in 
the analysis of the optical spectra) might account for a factor of 
three in the difference only. A combination of several other factors 
can explain the rest of the difference : (i) In these weak wind stars, 
the influence of porosity/vorosity effects, and the possible presence 
of hot gas, can potentially bias the estimated mass-loss rate by fac¬ 
tors from 3 to 10; (ii) In the region of parameter space correspond¬ 
ing to these stars, CMFGEN and FASTWIND produce system¬ 
atically different results for optical lines, with more emission/infill 
predicted by CMFGEN. 

• For the re-visited metallicity, we find that the mass-loss rates 
of IC1613-A13, the earliest star of the sample, is in good agreement 
with the mass-loss predicted with the recipe by | Vink et al.| ( |20OT| l. 
The C IV P Cygni profile of A13, which is stronger than that of the 
SMC dwarf MPG 324 (cf. Fig.[^ is probably a consequence of its 
more evolved status. 

• On the other hand, M measured with FUV spectra for IC1613- 
B11 and WLM-Al 1, two late-type O supergiants, are clearly lower 
than measured in the optical by |Tramper et al.| ( |2014^ . Because of 
the systematic errors, they may however be consistent with the the¬ 
oretical predictions by |Vink et ^ j2001| > for the preferred metal- 
licities. If the later is true, porosity/vorosity effects and the state of 
the gas (i.e., what fraction is hot) are of crucial importance. 

• We used a different approach to calculate mass-loss rates, as 
outlined recently in |Lucy| ( |2012| l. With this alternative method, 
avoiding scaling laws and assumptions on the terminal wind ve¬ 
locities, the measured mass-loss rates of the three stars are then in 
much better agreement with theoretical values. We therefore con¬ 
clude that there is no indication at this moment that the radiatively- 
driven wind theory breaks down at metallicities at or slightly less 
than SMC values. 

Our results are restricted to a very limed number of objects. 
We pointed out that the observed mass-loss rates within a spec¬ 
tral/luminosity class may vary within a factor of four or more just 
because of the scatter of the intrinsic individual parameters. To 
reach a firm conclusion on the nature of the M(Z) relation there¬ 
fore requires larger statistical sample, preferentially observed in 
the FUV at (relatively) high spectral resolution and high signal- 
to-noise ratio, to simultaneously derive the iron content and the 
wind parameters. This, unfortunately, is beyond routine capability 
of HST for galaxies more distant than the MCs. On the other hand, 
a fairly decent sample of objects is available in the HST archives 
for stars in the SMC and we shall now focus on these objects to 
study the properties of stellar winds at low metallicity. We engaged 
in such a study and results will be presented in a coming paper 


discussing the properties of mass-loss for stars from dwarfs to su¬ 
pergiants. 
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APPENDIX A: BEST FITS 

In this appendix we compare our best-fit models to the optical 
and UV spectra. Absolute fluxes in UV spectra are provided af¬ 
ter correcting for reddening (see Sect. ® and are expressed in 

A — 1 —1 —2 

. s . cm 
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Wavelength (A) 


Figure Al. FUV spectrum of IC 1613-A13 (in black). Oveiplotted in red is the best fit model (see text for comments). A satisfactory fit is obtained for the 
observed photospheric (mainly Fe forest) and wind profiles. In particular, the N V A1240 and N IV A1718 wind profiles are well reproduced, indicating con'ect 
ionization factions in the wind. 
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Figure A2. Optical spectrum of IC 1613-A13. Overplotted in red is the best fit model (see text for comments). 
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Figure A3. FUV spectrum of IC 1613-B11 (in black). Overplotted in red is the best fit model (see text for comments). 
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Figure A4. Optical spectrum of IC 1613-Bl 1. Overplotted in red is the best fit model (see text for comments). 
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Figure A5. FUV spectrum of IC WLM-Al 1 (in black). Overplotted in red is the best fit model (see text for comments). 


© 2013 RAS, MNRAS 000.[T](^ 










































Normalized Flux 


26 


J.-C. Bouret et al 



3900 4000 4100 4200 4300 4400 



4830 


4855 


4740 


5100 


4875 


6510 6528 6546 6564 


6582 6600 


Wavelength (A) 


Figure A6. Optical spectrum of WLM-Al 1. Overplotted in red is the best fit model (see text for comments). 
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